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ABS TRACT 


This thesis describes the design and construction of a wide band 
(60-3060 Hz), three phase (including zero sequence), non-transposed, 
lumped element, power transmission line model. The model consists of 
from 15 to 65 “three phase L sections" - depending upon the length of 
the actual line being modeled - and includes trimming features to model 
both horizontal and vertical configurations of lines in the 345 kV to 
1500 kV range. 


To predict the performance of the model and determine the number 
of sections required to model a given line, the poles of the driving 
point impedance of various numbers of L-sections in series is determined 
by computer and compared with the corresponding poles of the actual line 
being modeled. For either an open or short circuited termination, the 
performance specification is that in the range from 60 to 400 miles long, 
the percent error between the poles of the model and those of the actual 
line shall not exceed 5 percent. This criterion results in the require- 
ment for 15 to 65 L-sections. Thirty sections are of 4 mile "length" 
and thirty five are of 8 mile "length". 


With the number and "length" of the model sections determined, 
electrical parameter data for a number of typical extra high voltage 
(EHV) lines are suitably scaled to select a base value of model inductance 
and capacitance for the 4 mile and 8 mile sections. The effects of non- 
transposition are modeled by incorporating two mutual inductors between 
the phases. Zero-sequence inductance is modeled by a specially designed 
low Q inductor whose inductance decreases appropriately with increasing 
frequency. 


To ensure that the overall quality factor of the model remains 
greater than or equal to that of actual lines for all frequencies between 
60 Hz and 3060 Hz; low loss, powdered molybdenum-nickel-iron toroid cores 
are selected for tne inductors, and low dissipation factor polystyrene 
capacitors form the basic capacitance of each section. With the designed 
quality factor of the model exceeding that of most EHV lines in operation, 
the Q of the model is brought into correspondence with actual lines (at 
60 Hz) by suitably inserting resistance into each section. 
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CHAPTER 1 
INTRODUCTION 
1.1 Transmission System Simulators 
The three phase power transmission line model to be designed 
will form part of a power transmission system simulator (TSS) for use 
in the W.1.T. Electric Power Systems Engineering Laboratory. Knowledge 
of the transient behavior of power system networks is essential to the 
system designer and planner in order to correctly dimension or rate 
the various components of a power system. Circuit breakers and 
insulators, for example, cannot be correctly sized without some 
knowledge of the magnitude and duration of system fault currents or 
Switching overvoltages. With 765 kV systems already in operation, 
and 1100 kV and 1500 kV systems proposed for the future, it is 
particularly important to know the magnitude and shape of over- 
voltage transients caused by routine switching operations. Models 
for 1100-1500 kV lines still in the design stage are of obvious value, 
and the very limited access to extra high voltage (EHV) lines in 
operation makes a model of these lines highly desirable. 
1.2 General Model Parameters and Constraints 
A transmission line model (t.1.m.) is only one part of a TSS. 
Most simulators include, in addition, models of power system generators, 
transformers, compensating reactors, lightning arresters, circuit 
breakers, ground fault neutralizers, etc. In the present t.I.m. 
design, the requirement for interface with other TSS components places 


certain constraints upon the model. In particular, compatibility 








les 


with existing TSS transformers dictates that the nominal operating 
voltage of the t.].m. be 138V line-to-line and that nominal phase 
current be 38.6 ma. Other required features for the t.].m. are that 
it be capable of modeling vertical, horizontal, transposed, and non- 
transposed EHV line configurations from 60 to 400 miles in length, 
have quality factors exceeding those of existing EHV lines from 60 
to 3060 Hz, model the variation of zero sequence inductance with 
frequency over the range 60 to 3060 Hz, and model lines in the 
voltage range from 345 kV to 1500 kV. In addition, the positive 
Sequence inductance and capacitance of the t.].m. must remain 
essentially constant over the range 60 Hz to 3060 Hz, and inductance 
must remain constant with varying current in the range from 0 to at 
least one ampere-thus ensuring constant model inductance up to 
25 p.u. fault current. 
1.3 Features and Scope of Design 

One of the more unique features of this t.].m. design is 
that it is to be capable of modeling non-transposed power 
transmission lines. Numerous transmission line models in simulators 
are presently in operation! ~*, but all in the U.S., at least, are 
transposed models. A literature search covering the period subsequent 
to 1960 disclosed only a limited literature on t.1].m.'‘s - all 
relating to transposed models. Many actual lines are, in fact, not 
transposed. Previously, with relatively lower voltage transmission 
lines, it was assumed that the small effects of non transposition 
could be ignored. Non-transposed lines were mathematically and 


physically modeled as though they were transposed, and model results 





were considered to agree satisfactorily with actual lines. Recent 


9 on EHV lines indicate, however, that the effects of non- 


studies 
transposition can no longer be ignored. This thesis is an attempt 
to design a t.].m. which will include the effects of non-transposition 
and document the design in sufficient detail that persons other than 
the designer can construct and operate the model. 
1.4 Extensions of the Design Concept 

The t.l.m. design discussed in this thesis is intended 
specifically for the modeling of electric utility power transmission 
systems for which the X of Xp ratio (per unit length) at 60 Hz is of 
the order of 230,000-2/0,000 mic (X 6 = positive-sequence capacitive 
reactance; Xe = positive-sequence inductive reactance). The design 
concept, however, 1S not at all limited to utility lines. As will 
become clear in later sections of the discussion, any band of 
Xo/X ratios could have been provided - had it been desirable or 
necessary. With the proper selection of model components, any audio 
range power transmission system involving inductive and capacitive 
coupling could be accommodated. Thus, the model concept is equally 
Suited to modeling, for example, under ground power cables, long 
ship-board or industrial plant cable runs or long synchro cable runs, 
etc.. In addition, the mathematical approach developed in Chapter 3 
for the solution of the impedance polynomial has application in any 


area where the impedance concept has validity - including, for 


example, the mechanical impedance in mechanical vibrations problems. 
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CHAPTER 2 
MODEL CONFIGURATIONS 
2.1 General Non-Transposed Models 
A section of a general, lJumped-element, three-phase 
transmission line is Shown in Figure 2.1]. The electrical para- 
meters of the model follow directly from the inductance, 
capacitance, and resistance matrices of the actual line. With the 
liatrices of the line defined as 


Laa Lab Lac 


LL] = Lab Lob Lic : 


Lac Lic bic 


(2-1) 
Caa “Cab “Coe 
Be ere sy seme 
“C56 “Che Coe 
Raa Rab Rac 
and [R] = Sols Rib Ric ; 
R R R 


dG = bCaw ec 
tne correspondence between line and model inductance and capacitance 


paralieters is as follows (See Appendix A.1): 


Loa 7 =i) 7 LOO Lab Li2 " Loo 

Lyb * 422 * Lao Lbe * 423 * Yop 

aga 8 © Loe eae. Iaeees00 (2-2) 
bap = “12 So on ae lone 

Femur Gob | 12 e238)” “20 

re ~ 693 Coe ~ C13 b53 i C39 
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Tne resistance Rog in the model represents only an approximate model- 
ing of line inter-phase resistances. Rather than providing three 
separate inter-phase resistance circuits in the model, each of the 
three line inter-phase resistances is approximated by the average 


of the three ij.e., 


iM v ac _ : 
Rab ~ Bho = Rac g aa Ry (2-3) 
So that, 
Raa on 
AV, 
LR] = Ry Rob Rm (2-4) 
Ri Rn Ree 


Comparing line equations witn those of Appendix (A.1), it is found 


that: 


Rog = 8, 


Raa” Ria * ce 
Rob ~ Roa + Ry 
Roe : R33 i Rn 


Tne line paraiieters given in equation (2-1) include the effect 
of ground wires, so that the model parameters also include ground 
wire effects. The mutual inductances shown in Figure 2.1 represent 
“ideal mutual inductors adding no self inductance to the lines. 
The value of the mutual inductances is nominally about ten percent 
or less of tne self inductances. 

The circuit of Figure 2.1 can be transformed by simple 
algebra (See Appendix A.2) to the configuration of Figure 2.2. 
The circuit of Figure 2.2 is as general as that of Figure 2.1] but 


has one less mutual inductor. In Figure 2.2 the transformed 
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quantities, in tenns of Figure 2.1, are: 


oes boa = L ee tBu 


ais Cin 473 12 7 442 - '43 10 = “10 
Log = boo - by3 Log = bo3 - by3 C50 = C29 
Lay = yo Loo = Loo * £43 C30 = %30 
C2 = %y2 RS) en Roo = Roo 
ae 513 Rao = Roo 
53 = G93 Nom hos (A~4) 


If the line being modeled is of horizontal configuration; then for 


both circuits: 


Lig = lag - by) = bog = (lag - by3) 

bay = (Lay - by3) = b33 = (b33 - by3) 
10 “10 "30 + 30 

Cia = Cin = C53 = O53 

‘Twice ation = ge 


The model section shown in Figure 2.2 (hereafter referred to as 

“the model") is the configuration selected for the design. While 

it is by no means the only configuration possible, this model has 

a nuinber of advantages over otner possible inodels of equal generality. 
First, it should be pointed out that in order to economize on the 
number of cores to be purchased, the mutual inductors are wound over 
tne self inductors in the line. The configuration of Figure 2.2 
requires four cores, and incorporates the feature that no Self- 
inductor will be required to physically accept more than one mutual 
winding (the self winding leaves sufficient window area for only one 
Iutual). Second, the configuration in Figure 2.2 ensures that no 


phase line is required to incorporate the resistance of more than 
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one mutual winding. (The resistance of each phase line must be 
kept to a iinimum to keep Q high.) Third, when modeling horizontal 
lines, the mutual inductances of the Figure 2.2 configuration are 
equal-thus incorporating into the model the symmetry which actually 
exists in the line being modeled. (This would not be true of some 


other non-transposed models.°) 


2.2 Transposed Models 

As wentioned in Chapter 1, it is common practice to simplify 
the analysis of non-transposed lines by assuming that they are 
transposed. In most non-transposed lines the values of the three self 
inductances of the phase conductors are approximately the same, and the 
values of all inter-pnase parameters are not too different from one 
another. As shown in Appendix A.3, the transformation to a transposed 
View-point "averages" the line parameters and "re-distributes" them 
uniformly among the three phases. The result is a balanced, symmetrical 
model which for many applications is a good approximation to what is 
in reality a non-transposed line. A schematic of a transposed line 
model section resulting from the transformation of Appendix A.3 is 
Shown in Figure 2.3. In terms of the general model of Figure 2.1, 


the parameters of Figure 2.3 are: 


a = “10 * “20 * 30 (A-5c) 
ah z “ia * 23 * "13 (A-5d) 
Rune "u* Wee * O33 (A-5e) 
e = ru * ee * ©33 (A-5a) 
e = nia * £23 * 713 : (A-5b ) 
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Zc0el Symmetrical Component Equivalent Circuits for Transposed Model 


In Appendix A.3 1t 1s shown that the transposed line model 
can be reduced to two symmetrical component models. The results are 


repeated here for continuity. With the definitions 


oT 7 oe 
Ros, = Rs RI og, = Rg + SRoq pig, = Cg + Sly 
Ff ST 
ees. 7 <0 Lp sg, F bo - by 
eae 
and Las. = Lo + él. + SLi ys 


the transposed line equations (A-7) - with positive sequence currents 


only-reduce to: . 
- is R I +t aa 
OZ Pesce Pes oe 
- a R IL +k *’b (A-8) 
dZ Ploae DD Pyoee Oe 
OV. a1. 
- ~=—=R mae © 


RY PS 3c PS. ot 
And The Current Law equations (A-11) reduce to: 


dl oV 


e “a 
wove = Cbs, ot 
9] Vv 
b _ b 
Se Sas (A-12) 
zs ee = c eel 
OZ PeStrot. s 


These results show that for positive (or negative) sequence voltages 
and currents, the transposed line may be viewed as three identical 
Single phase circuits (the 120° phase difference understood). 
Similar conclusions apply to zero sequence voltages and currents 
where it can be shown (Appendix A.3) that under these conditions 


the transposed line equations reduce to: 
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a. ml 
epg) Ge a SF 
3V aI 
b b 
ape 2 Nes ye ey) Se (A-9) 
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oe Suse be us Loos ae 
and 
Dee A 
OZ Z.S. at 
aI 3V 
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pen? Ge (A-13) 
; 3V 
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Thus, the non-transposed model can bé approximated by a transposed 
model, and the transposed model can, in turn, be analyzed into two 
Symmetrical component equivalent circuits - one for positive-sequence 


and another for zero-sequence. 
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CHAPTER 3 
NATURAL FREQUENCY ANALYSIS 


3.1 Methods and Criteria for Predicting Model and Actual 
Line Correspondence 


As discussed in Cnapter 2 and Appendix A, the transposed 
line can be analyzed into two, symmetrical component, single-phase 
equivalent circuits. A single-phase L-section model as "seen" by 
positive (or negative) sequence voltages and currents is shown in 
Figure 3.la. A similar L-section model as seen by zero-sequence 
quantities is shown in Figure 3.1]b. Figures 3.la and 3.1b reflect 
the basic "L" configuration ifiptied by the tnree-phase models of 
Figures 2.1, 2.2 and 2.3, but "mn" or "T" section models are equally 
valid. In fact, under certain circumstances and with certain loads, 
the "nr" or "T" configurations may have advantages over the "L" 
models - these circumstances will be discussed more fully in later 
sections. 

One of the purposes for which the transposed line 
approxiniation can be used is that of establishing a "math-model- 
predictor" for conparison of an actual line and its proposed models. 
Some of the early questions in the design of a lumped element 
t.].m. are: How many sections are required? How "accurate" is the 
model with ten-twenty-forty sections? Should the sections be "L" 
or "nm"? What should be the "length" of each model section? For 
wnat frequency band must the model “correspond” to the actual line? 
Exactly what is it that should "correspond" between model and line- 
to what "accuracy"? A method of theoretically comparing an actual 


line and its proposed models is clearly desirable. If a t.].m. 
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could be constructed to exactly duplicate in magnitude and phase the 
impedance function of the actual line for all frequencies, then the 
iodel could be said to have "perfect accuracy". Perfect accuracy 
would probably require more lumped element sections than economically 
practical, so that some compromise is required. 

As tne first, and perhaps the simplest, approach toward 
answering the questions above, this design uses aS a comparison 
"device" the poles of the driving point impedance functions. The 
poles of a given length of actual line are compared to the poles of 
a proposed m-section wwodel. The number of sections, m, in the model 
is "adjusted" until tne frequency difference between corresponding 
poles of the niodel and actual line is sufficiently small. In each 
case, the actual line and model compared are the transposed, 
symmetrical-component equivalents of the three-phase systems. This 
simplification is sufficient for design analysis since it is not the 
detailed nature (e.g. “splitting") of the poles that is sought, but 
rather, a reasonably accurate estimate of their mean position in the 
frequency spectrum. Tie mean positions can then be compared to 
estimate the difference between the actual line and model. The 
transposed, synmetrical-component approximation is applied to both 
the matnematical model of the actual line and that of the t.l.m., 

So that errors introduced by this approximation are self compensating. 

In analyzing tne driving point impedance function of both 
the line and models, it is convenient and sufficient to consider 
two possible extremes as terminations. Both the sources and loads 
connected to the terminals of the line are non-linear devices (at 


least under some conditions of operation or fault) so that, as 








“ai 


viewed from the line, the terminals of the line may simultaneously 
appear as a high impedance to some transient components and a low 
impedance to other components. This "high-low" terminal impedance 
characteristic is provided for in the analysis to follow by consider- 
ing both open and short circuit terminations. The decision as to 
what is a "sufficiently" siall difference between the line and the 
model is deferred to a point later in the analysis where it can be 
lade to appear a little less arbitrary. 
3.2 Driving Point Impedance of Continuous Lines 

The relations describing the driving point impedance of 
continuous lines are well] known and will not be derived for this 
discussion. Readers interested in the details of the derivation 
inay consult any good text on transmission lines such as references 
7, 8, or 9. For the general transmission line circuit shown in 


Figure 3.2, the pertinent impedance relations for a length, &, are: 





Lec. = Z, tanh y2 (3-1) 
C= Z, coth yz. (3-2) 
wnere 
Lec. = Impedance at Z = & with a short 
circuit termination. 
Loc. = Impedance at Z = & with an open 
circuit termination 
Se at ae 


y= Y(R + sL)sc (3-4) 
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Ss = Complex piase constant 

R = Positive or zero-sequence resistance per unit length 
L = Positive or zero-sequence inductance per unit length 
C = Positive or zero-sequence capacitance per unit length 


Tne shunt conductance, G, nas been ignored in these relations. In 
tne short circuit case, tne poles are seen to occur when 
cosh y& = 0. (3-5) 
Setting y = ad jy; » equation (3-5) can be expanded as follows: 
cosh(y + Jy; )% = 0 
cosh(y .2)cosh( jy.) + sinh(y 2)sinh(jy;2) = 0 
cosh(y 2) cos YR +4 sinh(y2)sin Y;2 = 0. 


Thus, for cosh y& = 0 it is necessary that 


cosh(y 2) cos yi2 = 0, (3-6) 
and 

sinh(y2)sin yi2 = 0 (327) 
Case No. 1: 

Suppose ees Q 

Then eq. (3-7) is satisfied and eq. (3-6) 

implies y = jy; = + j(2ntl)1/2 noOe ls cee ee (3-8) 
Case No. 2: 

Suppose ven # 0 but y;2 = +> n=0,1,2.... 


Then eqs. (3-6) and (3-7) are never both satisfied. 
Therefore Case No. 1 gives the only solutions for y which result 
in the poles of equation (3-1). A similar development applied to 
equation (3-2) gives for the open line: 


yQ = + jnn eee eee (3-9) 
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If equations (3-8) and (3-9) are expanded and solved for the complex 


phase constant, s, tnere results: 


2 ame 
: ae oe +/ 6D SGV 2 (3-10) 


“Al 


R R C n° Tt 
oe lap! 2 Gap) = = n=0,1,2.... (ai) 





— 
W 


If the line shown in Figure 3.2 is lossloss (R = Q) then: 


Z. .. = JZ, tan Be (2212) 

Die = - jZ, cot BL (3-13) 
_ i 

Zz, =/¢ (3514) 

8 = wvLC = phase constant (3-15) 

Wee > eel aD n=0,1,2.... (3-16) 
eal 22 VLC 

W ol n=0,1,2.... (Bata 


O.C. me) JG 


3.3 Driving Point Impedance of L-Section Models 





A schematic diagram of an L-section model for nodeling the 
line of Figure 3.2 1S shown in Figure 3.3. Both figures refer to 
the symuietrical component equivalent circuits of their respective 
three-phase systems. Thus, both figures represent a "transposed- 
line-approximation" to their respective non-transposed, three-phase 
counterparts. In the case that the L-section line is shorted, a 
clear distinction must be made as to which end is the source and 
wnici} is the termination. In the first case, a short can be 
placed across tne last capacitor-removing it from the circuit. This 
case is referred to as “inductor input". In the second case, the 
Short can be placed to connect the "last" inductor to ground. This 


case 1S referred to as "capacitor input". 
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3.3.1 Inductor Input L-Section ilodel (Shorted) 


The inductor input case is that shown in Figure 3.3. The 
derivation of the driving point impedance expression is rather 
lengthy and is shown in detail in Appendix B.1]. The lumped element 
parameters of Figure 3.3 are defined in terms of the distributed 


parameters of Figure 3.2 as: 


REI 

R. a (3-18) 
Sree 

Le =i (3-19) 
muy 

C. Shei (3-20) 


ye | 


Where "m" is the number of sections in the model. The result from 
Appendix B.1 for an m-section; lossy model is: 


(R.+L, )sinh mA 


‘sec. ~ sinh mi-Sinh(mT)d * (B-16 ) 
The denominator of equation (B-16) is a polynomial in the real 


variable cosh A so that the solutions, [cosh A] to the equation 


POCO NE 
sinh mA-sinh(m-1)\ = 0 yield the propagation constants in the 


propagation relation 


ys SS 


[cosh Meee = 5 aaa (B-7) 








Solving equation (B-7) for s gives: 


/ @ 2 
-sic. Gp us Go) - “Z-< [1-(cosh More (B- 26) 


If the model of Figure 3.3 is lossless (R. = 0), 
Jul. Sinh ma 


“sic. ~ Sinn mAsstnan(m=T)A MM V+223---. (B-17) 
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wb C, 
cosh A = 1 - 5 (B-18) 
ee eet A), tic) 2.2 (B-27) 
Shek Q LC root : lg 


For an example application, suppose their are seven sections 
in a proposed model. Then m=7; there are six values of [cosh Ras 


and six values of s or 


S.C. Ooo The six poles of the model given by 


equations (B-26) or (B-27) are compared to the first six poles of the 
actual line given by equations (3-10) or (3-16). The difference in 
frequency between two corresponding poles is a measure of the error 
of the model for that pole number. Equations (B-16), (B-26), and 
(3-10) were implemented on the computer, as shown in Appendices C and 
D, for various numbers of sections up to about forty. A plot of the 
resulting pole error vs. pole number is shown in Figure 3.4a. The 
percent error in Figure 3.4a is defined as 
[ equation (B-26)-I [equation (3-10) ] 

% Error = 100 ~——T [equation (3-TO)] : (3-21) 
The parameters used in equations (B-26) and (3-10) for Figure 3.4a 
are the positive sequence values of a 765 kV “typical aluminum 


tower" configuration: 


2 = 100 miles (3-22) 
R = .0251 ohms/mi (3-23) 
C = .021 uf/mi (3-24) 
L = 1.45 mh/mi . (3-25) 


The pole number shown in Figure 3.4a is, of course, a discreet integer 
variable with pole no. 1 corresponding to the pole of lowest 


frequency in either the actual line or model. The curve for a given 
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Plot of Pole Error vs. Pole io. 
L-Section Line (Pos. Seq.) 
Inductor [Input 


: Shorted Termination 
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21 Sections 


—— 39 Sections 


Figure 3.4a 
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number of sections in the model is a series of discreet points but is 
drawn continuous to avoid confusion. Figure 3.4a shows, for example, 
that at least twelve sections are required to model the line if it is 
desired to keep all of the first five poles of the model within five 
percent error. 

If the typical parameter values of equations (3-22)-(3-25) 
are substituted directly into equations (3-10) or (B-26), it is found 


that the term > is of order 10' whereas the terms 





2 2 2 
(2nt+1)o1 2m ’ 
a or cd [1-cosh eer 
7 


are of order 10° or 10° for even the lowest values of moron. This 
Suggests that for purposes of a pole error analysis, the resistance 
terms can be neglected when considering the positive sequence of the 
lines. This, in turn suggests that a pole error analysis using the 
lossless equations (B-17), (B-27), and (3-16) should give results 
very similar to the lossy case shown in Figure 3.4a. The results 


for the lossless case paralleling the analysis leading to Figure 


3.4a are shown in Figure 3.4b. The percent error in Figure 3.4b 


es [1-(cosh i) i (én+1)q 
R Lie root oe VLC 


% Error = 100 — (entl)m OSS (3-26 ) 
22 vLC 


is defined as 


As expected, Figures 3.4a and 3.4b are practically identical so that 
losses will be ignored throughout the remainder of the pole error 
analysis. flote that equation (3-26) is independent of 2, L, and C, 
SO that the pole error curves of Figure 3.4b apply to any line and 


model provided the model is “inductor input", L-section, and both 
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Plot of Pole Error vs. Pole No. 
L-Section Line (Pos. Sea.) 
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Figure 3.4b 
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line and model are terminated in a short. The pole frequencies 

corresponding to the first fifteen poles of a shorted, typical 

765 kV line, for various line lengths are shown in Table 3.1. 
Equations (B-26) and (3-10) for the configuration of Figure 

3.3 were also run uSing zero sequence values of a typical 100 mile 

765 kV line. The pole error results for this case are shown in 


Figure 3.5 with 


& = 100 miles (3-27 
R = .51 ohms/mile (3-28 
C = .0143 uf/mile (3-29 
L = 4.2 mh/mile . (3-30 


The percent error in Figure 3.5 is as defined by equation (3-21). A 
comparison of Figures 3.4a and 3.5 shows that the pole error for zero 
Sequence parameters is practically identical to the error for positive 
sequence parameters. In addition, Figure 3.5 is practically identical 
with the lossless results shown in Figure 3.4b. Thus, throughout the 
remainder of the pole error analysis, the lossy line and zero sequence 
pole error will be assumed to be the same as the lossless line error 
of any given configuration. 

3.3.2 Capacitor Input L-Section Model (Shorted) 

If the model of Figure 3.3 is turned around so that the 
Capacitor is at the input, the result (lossless) is shown in Figure 
3.6. The derivation of the driving point impedance of this line is 
carried out in Appendix B.2. The analysis considers only the lossless 
case since, as discussed in section 3.3.1, the effect of resistance 


can be neglected in the pole error analysis. The results from 
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Table 3.1 


Short Circuit Pole Frequencies (Pos. Seq.) 


Pole 

No. 60 Mi 100 Mi 200 Mi 300 Mi 400 !1i 
] 760 Hz 456 Hz 228 Hz LS2 Hz 114 Hz 
2 Zecou lies70 685 456 343 
3 3,800 2,280 1,140 760 570 
4 5,345 3,200 1,600 1,054 800 
5 6,840 4,100 2,050 1,370 15025 
6 S905 5 ,020 2,560 1,675 13258 

7 9,850 5,920 2,960 i975 1,480 
8 ie oo0 6,840 3 5420 Z5Gh9 Varo 
9 12,900 75750 3 ,880 2,580 1,940 
10 14,400 8,650 4,325 2,880 2,160 
1] 15,950 9,580 4,790 3,190 2,390 
12 17,410 10,470 asks 3,490 2,620 
13 18,950 11,390 5,695 3,790 2,840 
14 20 ,500 12,300 6,150 4,100 33075 
15 22,000 13,200 6 ,600 4,400 3,300 


Table of Pole Frequencies for Shorted, Typical 


765 kV Line (Positive Sequence-Lossless) 
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Plot of Pole Error vs. Pole No. 
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Figure 3.5 
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Appendix B.2 are, 


Jul. Sinh mA 


25.c. ~ Sinh(mT)a-sinn my MT 9ee oe sedis 
wl C, 
cosh A= 1 - =) a (B-33) 
Ww = a [1-(cosh A) ] m=1,2,3 (B-40) 
CC amy LG Root Sire eae 


Note that equation (B-39) is not the same as equation (B-17), so that 


the values of (cosh A) obtained from (B-39) for use in (B-40) will 


root 
be different from the roots obtained from (B-17) for use in (B-27). 
Equations (B-39), (B-40), and (3-16) were solved with the computer 

as shown in Appendices C and D. A plot of the resulting pole error 
vs. pole number 1s shown in Figure 3.7. Having neglected resistance, 
the poles of the model are given by equation (B-40) and those of 

the actual line by equation (3-16). The percent error in Figure 3.7 
is defined by equation (3-26). Figure 3.7 shows that the error 
behavior for the capacitor input case is radically different from 

that of the inductor input. The capacitor input model consistently 
under-estimates the actual pole frequency, whereas the inductor input 
model over-estimates the frequency of the lower order poles, and, 
throughout its range, yields pole frequencies higher than the 
Capacitor input model. The inductor input model gives higher pole 
frequencies due to the fact that the shorted termination removes 

the last capacitor from the model-thus reducing overall] model capacity 
and increasing the model's natural frequencies. Note also, that the 
Capacitor input model errors are larger at lower order poles and 


smaller at higher order poles than corresponding lengths of inductor 
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input model. 
3.3.3 L-Section Model, Open Termination 

In the case of an open termination, the poles are unaffected 
by which terminals are considered to be the input or output, and only 
One analysis is required. The schematic of the open circuit case is 
just Figure 3.3 with the short removed. The driving point impedance 


for this case is derived in Appendix B.3. The results are (R. = 0): 


7 Jeet sinh(m+1)A-sinh mA 
o.c. Ys 12 Sinh ma-sinh(m+])A-sinh(m-1)a (B-46) 
wel C. 
cosh A = 1 - aa (B-18) 
2M [2 7y_ b 
yc.  2./ LC [1-(cosh A) J (B-47) 
Where, again, use has been made of the fact that, 
Slee Ce 

ts am? “Ss me Ss 

m = number of sections in the model, 

£ = length of the line being modeled, 


and L and C are the electrical 

parameters per unit length of any 

given line. 
Equations (B-46), (B-47), and (3-17) were solved with the aid of a 
computer as discussed in Appendices C and D. The results are 


shown in Figure 3.8 with the pole error defined as 





m 2 nt 
& fT LI-(cosh Moot - 7 
% Error = 100 a (3-31) 
gvLCc 


As in equation (3-26), (3-31) is independent of 2, L, and C, and 
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Plot of Pole Error vs. Pole No. 
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Figure 3.8 
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depends only upon line configuration. Again, the values of the roots 
are different from what they were for either of the preceeding cases. 
The pole frequencies corresponding to the first fifteen poles of an 
open, typical 7/65 kV line, for various line lengths is shown in 
Table 3.2. 
3.4 Driving Point Impedance of m-Section Models 

Referring to the line shown in Figure 3.6, it is seen that 
the L-section model can be converted to a 7 model by a simple 
adjustment of the first and last sections. If the capacitor at 
the mth node is replaced by one of value C/2, and another capacitor 
of value C./2 is added at the "O"th node, the result is a m-section 
model. A similar type of adjustment could be made using inductors 
of value L./2, and the result would be a T-section model. The t 
and T-models are reciprocal networks, whereas the L-model is not. 
A pole error analysis for m-sections follows in sections 3.4.1 and 
3.4.2. The analysis for T-sections is not included since it is 
assumed that it will add little to all the cases considered. 
3.4.1 m-Section Model (Shorted 

The derivation of the driving point impedance relations for 
the shorted, m-model are carried out in Appendix B.4 and the solution 
is carried out in Appendices C and D. The result for the shorted 
t-model is shown in Figure 3.9. The percent error in Figure 3.9 is 
defined as in equation (3-26), but again it should be noted that 


the values of [cosh i) differ for each line configuration and 


root 
termination. 





Pole 
NO. 


1] 
14 
13 
14 


60 Mi 
0 Hz 
a alte 
3,020 
4,540 
6 ,040 
75950 
9055 
10,580 
12,090 
13,600 
hS050 
16,600 
18,200 
19,650 
ZU 


Table of Pole Frequencies for Open, Typical 
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Table 3.2 


Open Circuit Pole Frequencies 


100 Mi 
0 Hz 


905 
1,810 
2,/20 
3,620 
4,540 
9 5440 
Ono oU 
fycou 
Sig JAS 
9,050 
Sieh: 


10,900 
11,800 
12,700 


200 MMi 
QO Hz 


452 

905 
1,360 
1,810 
Cac) 
Cw) 
Sig Gis) 
Sine) 
4,092 
sey 40) 
4,999 
9 5450 
5 5900 
6,350 


300 fi 
QO Hz 


cane, 

604 

906 
lpyaus 
Tote 
Pcie 
Zl al 
2,420 
120 
3,070 
S,025 
3,640 
3,940 
4,240 


765 kV Line (Positive Sequence-Lossless) 


400 Mi 
0 Hz 


226 

402 

680 

oT Oe) 
5) 
Ieee. 
Reiss) 
1,815 
2,046 
2,260 
2,500 
Calico 
Zo Ow 
oe 
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Plot of Pole Error vs. Pole No. 
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3.4.2 m-Section Model (Open) 


The open circuit case is derived in Appendix B.5 and solved 
in Appendices C and D. The results are as shown in Figure 3.10. For 
Figure 3.10, the percent error is as defined in equation (3-31) with 


due recognition that the values of (cosh i) in this case are 


root 
different from their values in the opened, L-section line. 
3.5 Conclusions 

1. The results shown in Figures 3.4a, 3.4b, and 3.5 indicate 
that for typical EHV lines and models, the percent error between lines 
and models is largely independent of losses in either. Thus, the 
percent error for either positive or zero-sequence is the same for 
a given number of model sections and a given type of termination. 

2. The results indicated by Figures 3.4b and 3.7 show that 
due to the non-reciprocal character of the L-section model, the error 
between model and line depends greatly upon which terminals of the 
model are considered input or output. In fact, the capacitor input 
case is seen to be considerably inferior to the inductor input case 
for low order poles. 

3. Comparison of the shorted L-model with the shorted 1- 
model shows that the L-model has a rather more complicated error 
pattern than the m-model. The L-model first over estimates the 
true pole position for low order poles - then under estimates for 
higher order. The n-model consistently under estimates the position 
of the actual pole. 

4. The open circuited m and L-models give practically 
identical pole error plots. (Actually, each m-model has one 


additional pole higher than the highest pole of the corresponding 
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Plot of Pole Error vs. Pole No. 
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L-model. ) 

5. Any L-section model can be converted to either a 7m or T 
model by a simple change in the first and last sections to obtain 
any of the pole-error patterns of preceeding sections. 

6. A suggested method of using the pole error plots is 
Shown in Figure 3.11. Figure 3.11 shows, as an example, the pole 
error “curves' of a twelve section L-model as taken from Figures 3.4b 
and 3.8. A five percent error bound is suggested by the dotted lines 
at the top and bottom of the figure. The region between the solid 
curves is the “error band" of the twelve section L-model. If the 
terminals of the line appear as a predominantly high impedance to 
frequencies corresponding to the first five poles, the line error 
is given (approximately) by the "open" curve. A corresponding 
observation applies for the low impedance case. The dashed curve 
mid-way between the solid curves is meant to indicate the “nominal 
error’ when the terminals appear as neither open nor shorted. If 
twelve sections were used to model a sixty mile line, Figure 3.1] 
shows that the error band lies within the five percent error 
bounds up to about the fifth pole. As seen from Tables 3.1 and 3.2 
the fifth pole corresponds to about 6000 Hz for a sixty mile line. 
One could say then - “the 12-section L-model is good to within five 
percent up to 6000 Hz" (for a 765 kV line). An error plot 
Similar to that of Figure 3.11 could be made for any model type 
composed of any number of sections. In particular, if it is 
desired to model a 400 mile line to “within five percent" an 


inspection of Figures 3.4b, and 3.8 shows that a model of about 
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forty sections will ensure that the line is “good to within five 
percent’ up to the fifteenth pole. Tables 3.1] and 3.2 show that 
the fifteenth pole corresponds to a frequency of about 3100 Hz 
(for a typical 765 kV line). Thus, a minimum of forty ten-mile 
sections would be required to model a 400 mile line to within the 
five percent criteria. 

Tables 3.1] and 3.2 give the pole frequencies, specifically, 
for a typical 765 kV line. Actually, as will be shown in Chapter 4, 
corresponding tables for most other EHV lines would be very similar 
to Tables 3.1 and 3.2 (within about ten percent). With Tables 3.1 
and 3.2 as representative of EHV lines, it is now possible to make 
some quantitative decisions concerning the number of sections to be 
built for the t.].m. 

7. As stated earlier in this report - as a minimum - the t.1].m. 
is to mode] EHV lines from sixty to 400 miles in length. The five 
percent “specification” referred to earlier is actually a result of 
observations such as those in Conclusion 6. It was noted earlier 
that forty sections 1s the minimum number which will just ensure 
that a 400 mile line-model is within five percent to 3060 Hz. This 
in turn, implies at least forty ten-mile sections for the model. In 
addition, it is proposed to model all line lengths between sixty 
and 400 miles by changing the number of sections in the model. Thus, 
if a line to be modeled is 396 miles long, it would be desirable to 
set up a model "396 miles long". This adjustment could be 
accomplished with the ten-mile sections by using "taps" to adjust 


length. However, this t.].m. design must model, horizontal or 
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vertical, non-transposed lines - this too requires taps. To avoid 
having to tap the sections for length, it is proposed to make the 
sections "fine" enough (i.e. short enough) so as to make tapping 
unnecessary for length adjustments. Line lengths to 400 miles can 
be reached to within one percent by making the largest sections 
have only eight-mile length. Thus, the 396 mile line, e.g., could 
be reached (to one percent) by either forty nine or fifty eight- 
mile sections. But eight-mile sections are not "fine" enough to 
model short lines, etc., etc. Finally then, with considerations 
Such as those above, it was decided that the t.].m. design should 
consist of thirty four-mile sections and thirty-five eight-mile 
sections. This selection is seen to be quite conservative (by 
the chosen criteria) and provides a model somewhat "finer" than 
some other t.1.m.'s./74 

8. The pole error analysis is seen to provide a means of 
quantitatively comparing proposed models, and is useful in indicat- 
ing the relative merit of, say, doubling the number of sections - 
or halving them, etc. 
3.6 Applications and Recommendations 

The solution technique employed to solve the impedance 
functions of this chapter can have application in other fields of 
Study. The denominator of each of the impedance functions of this 
chapter corresponds to the "frequency polynomial" of mechanical 
vibrations problems. If the frequency polynomial is set equal to 
zero, it becomes the "frequency equation" for determining the 


natural frequencies of mechanical systems. Even in systems with 
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only two degrees of freedom the frequency equation is of fourth degree 
(quadratic in a) so that, in continuous systems or systems having 

more than two degrees of freedom, complicated solution techniques 

must be employed or approximations made. The technique suggested 

in this chapter and employed in Appendices C and D is a means of 

having the computer generate the frequency equation and then solve 

it for the natural frequencies of the system. A detailed investiga- 
tion into the exact class of problems for which this solution technique 
may have validity has not been conducted by the author, but it is 


suggested that such an investigation could be worthwhile. 
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CHAPTER 4 
PARAMETERS OF EHV LINES AND MODELS 

4.1 Parameters of Representative EHV Lines 

In order to determine the parameters or component values for 
a t.l.m., it is first necessary to know the parameters of typical 
lines the t.].m. would be expected to model. The t.].m. should be 
capable of modeling not only lines already in existence but, in 
addition, it should be capable of modeling proposed and even hypothetical 
lines. For information pertaining to the parameters of existing lines 
use was made of Report Wo. 710, Electrical Characteristics of Trans- 
mission Lines, of the American Electric Power Service Corporation 
(AEP) !°, Appendix IV of Report No. 710 describes in detail thirty- 
four different EHV transmission line configurations in use by the 
company. From the tnirty-four configurations a sample of $ix was 
selected as representative of types of configurations. For example, 
AEP type "H" was taken as one representative of a 345 kV, single 
conductor per phase, single ground wire, vertical configuration. 
Type "H-Special" was selected since it is a common line and represents 
a form of triangular configuration. A hypothetical type "39F" was 
"constructed" by using two phase-conductors per bundle spaced 
eighteen inches apart - this was taken as representative of a two- 
conductor, ‘vertical", 345 kV line. A type "41A" line was selected 
to represent a 345 kV, two-conductor, horizontal configuration. The 
AEP "Typical Aluminum Tower" configuration was taken as representative 


of most 765 kV lines (horizontal). A hypothetical 765 kV "H- 


Special" was constructed by scaling the 345 kV "H-Special" conductor 
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distances in the ratios that exist between the type "41A" and the 
"Typical Aluminum Tower". Finally, two hypothetical 1500 kV con- 
figurations were taken from reference 1] - a horizontal configuration 
(Geometry No. 2 and Bundle Drawing No. 4), and a triangular line 
(Geometry No. 3 and Bundle Drawing Ho. 4). 

Typical positive and zero-sequence parameters obtained for 
these representative lines are shown in Tables 4.1 and 4.2. The 
values shown in Tables 4.1 and 4.2 (except for Xo/X, and Q) were 
obtained with the aid of a computer using the “Electric Power Systems 
Engineering Laboratory (EPSEL) Line Parameters Program !¢ This 
program requires as input the geometrical configuration of the 
conductors of a given line (along with conductor and earth 
properties) and computes the inductance, capacitance, and resistance 
matrices of the actual line as a function of frequency. In addition, 
as well as the parameters of Tables 4.1 and 4.2, the program computes 
and prints the values of the model parameters of Figure 2.2. The 
model parameters are printed on the basis of a one-mile long section 
and assumes a scale factor of one between impedances of the line and 
model. Model parameter values at sixty Hz, on the basis described 
above, for the eight representative line models are shown in Table 
4.3 (earth resistivity is taken to be 100 ohms-meters). Referring 
to Table 4.2 note that, except for the 1500 kV lines, the positive- 
sequence (or average) Q at 60 Hz is about twenty or lower. The high 
Q values for the 1500 kV lines are questionable and would be 
difficult to duplicate in a practical model. Consequently, no 


attempt will be made to even approach the Q values suggested by the 
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1500 kV line values. With a Q of twenty, the 765 kV horizontal line 
will be taken as the minimum standard (in Q) to be achieved by this 
t.l.m. This means (referring to Table 4.2) that each of the sixty- 
five sections of the t.1].m. must have a positive sequence Q of at 
least twenty at sixty Hz. The positive and zero-sequence Q of 
transmission lines is not constant with frequency and, in fact, 
increases with frequency. Graphs of the positive and zero-sequence 
Q of representative lines are shown in Figures 4.] and 4.2. The 
positive-sequence Q is seen to increase steadily in the range to 
3060 Hz while the zero-sequence Q becomes nominally constant in the 
range from 1500-3060 Hz. Each model section must provide a Q 
characteristic equivalent to those of Figures 4.1 and 4.2 in order 
to accurately model the actual line. If the model section is 
designed to meet the positive-sequence Q requirement of the typical 
765 kV line, than it can be made to model lines of lower Q by 
simply inserting resistance in the phase conductors. 

The positive-sequence inductance of EHV transmission lines 
is practically constant over the frequency range of interest, but 
the zero-sequence inductance is not. A graph of the variation of 
zero sequence inductance with frequency is shown in Figure 4.3 for 
the typical 765 kV line. The zero-sequence inductance is seen to 
drop markedly with increasing frequency-particularly in the low 
frequency region. Thus, in the model also, the zero-sequence 
inductance must correctly decrease with increasing frequency and 


produce (as nearly as possible) the scaled equivalent of Figure 4.3. 
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4.2.1 Line/Model Scale Factors and General Model Parameters 

Due to the fact that the t.].m. must be compatible with existing 
TSS transformer and compensating reactor models, the base voltage of the 
t.].m. must be 138 V line-to-line, and the ratio of the base impedance 
of the model to that of a typical 765 kV line must be 0.355. This means 
that the base current and power of the t.].m. are 38.6 ma and 9.2 watts 
respectively. The base impedance of a typical /65 kV line on a 100 MVA 
base is 5850 ohms. The base impedance of the model is 2078 ohms. With 


the above impedance ratio constraint between the model and 765 kV line, 


le 

= model _ 55) 
B actual 

R 

~ model _ ace 
B actual 

and 

by mode |] ] 

CuneELNOGS «7 CoCo ° 
B actual i 


To model the typical 765 kV line, all the 765 kV parameters of Table 
4.3 must be appropriately scaled by the ratios above. For example, 
Loo of Table 4.3 becomes .564 mh/mi, or 2.26 mh/4 mi, or 4.52 mh/8 mi. 
Thus, the 4-mile section would be constructed with a basic 2.26 mh 
inductor. An identical procedure applied to the 345 kV type "41A" 


line gives (on a 100 VA base): 


(Z = 11902 


B) 345 


‘um, 2078 _ 1 44 
TZy)3q5 1190 


Loo = 3.1 mh/mi, or 12.4 mh/4 mi, etc. 
To model the 765 kV line the inductor in the four mile section must 
be 2.26 mh and to model the 345 kV line it must be 12.4 mh. This 


implies a very wide range of taps on the model inductors. A wide 
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range of taps on an inductor is not only difficult and expensive to 
construct, but has the further disadvantage that the Q of the inductor 
gets "tapped" as the inductance. Thus, tapping the inductance down to 
one-sixth of the total inductance cuts the Q down to about one-sixth 
of what it was for the total inductor. 

A procedure for solving this problem is to choose, for lines 
other than the typical 765 kV, a power base other than the commonly 
used 100 HVA. Again Suppose that the 345 kV type "41A" line is to be 
modeled. The first step is to establish an impedance ratio between 
the actual line and model such that Loo for the 345 kV model is the 
Same as that already provided for the typical (horizontal) 765 kV 


line. In other words, using Table 4.3, seek a scale factor F, 


pm (23) mode] 


(4-1) 
Zp line 
such that 
F(Lo0)3q5 = (Loo) 765 hor. * 2) 
So that 
L: 
22/765 hor. _ .564 
F = = = ,316, (4-3) 
boo) 345 1.78 
and 
Z 
pos 316 
Bg 
But 2 
1, = Be (4-4) 
BL PA 
so that 2 
(G3I6)aV,, 
pe eee SMA (4-5) 
BL Lee 


Thus, when modeling a 345 kV type "41A" line, it is necessary to use 
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a power base in the line of 18.1] HVA instead of the usual] 100 MVA, 
The resulting parameters required of a four-mile model section are 
Shown in Tables 4.4a and 4.4b respectively. Note in Table 4.4a that 
the inductors corresponding to phases "a" and "c" are equal - reflect- 
ing the symmetry of the horizontal lines. This is not the case for 
the "vertical" lines of Table 4.4b. With the inductances of the 765 
kV horizontal model as a base, note that all the L.'s of Table 4.4b 
are larger than 2.28 and all the L.'s are smaller. In fact, all the 
L,'s of Table 4.4b are larger than 2.28 by the same amount (A) that 
al] the L.'s of Table 4.4b are smaller than 2.28. The scale factor 
to achieve this "balance" for vertical configurations is not given 
by equation (4-3) but is determined as follows: 

Using Table 4.3, seek a scale factor F such that 


F(L = A 


Te vene. iz (L414) 765 hor. 


F ( - (L - A. 


L33) vert. A765 fa. 


Solving for F [using the fact that (L L 


11765 hor. “¢ 337765 hor. 2 


- vl 768 hor. si 
| 1] so tere 

The use of equation (4-6) ensures that the unbalance inherent in 

vertical lines is shared equally by both phases "a" and "b" in the 

model. A “balanced sharing" is desirable since a factor determined 

as in equation (4-3) gives the result that the phase "c" model 

inductor would be required to span a larger range than the phase "a" 

inductor - thus requiring more taps. Use of equation (4-6) ensures 


that the phase "a" and "c" inductors will have the same number of 


taps. Referring again to Table 4.4a, note that L for the two 
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Type 1500 kV 1500 kV 
AIA Hori zontal Triangular 
Zaey 2.34 2.16 
ence Z10 2.41 
2s 2.34 2.16 
ibe ke 0.630 eo 
0.27 0.23 0.11 
0.27 O23 0.11 
30.4 ZS «1 30 
30.4 28.1 30 
8.8 Ta 32 
166 164 12] 
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166 164 121 
65 10 1585 
73 6.5 19.8 
65 10 7.85 
230 67 70 


Table of Scaled Parameters Implied for a 
4-i4ile Model Section. (L in mh/4 mi., 
C in nf/4 mi., R in mQ2/4 mi.) 


Table 4.4a 
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Type 345 kV Type 765 kV 
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2.44 2500 Zoo 2.34 
2.34 Zell Yara) raya, 
Die ie 2.20 ngs re 
0.98 0.99 OR3 ele, 
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O20 -0.05 0.18 aN) Ole 
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Table of Scaled Parameters Implied for a 4-ilile 
Model Section. (L in mhn/4 mi., C in nf/4 mi., Rin m2/4 mi.) 


Table 4.4b 
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1500 kV lines is not 2.26 mh as would be found had equation (4-3) been 
used to determine the scale factor. In fact, still another technique 
has been employed to determine the scale factor of these two lines. 
This third technique will be discussed later in this chapter. 

Some remarks are in order at this point concerning the model 
zero-sequence inductance, Lis» and the inductor Ls of Table 4.4a. 

To obtain the zero-sequence characteristic required of a 4-mile model 
section for modeling a 765 kV line, for example, the graph of Figure 
4.3 is simply nied by the factor (0.355)(4.0) = 1.42; the result 

is shown in Figure 4.4. The model zero-sequence inductance is required 
to fall from 6.2 mh at 60 Hz to about 4.0 mh at 1000 Hz. The point 

to note is that L, of Tables 4.4a-b is not the zero-sequence inductance 
of the model per se, but L, 1s the component which must be designed 

to produce the characteristic shown in Figure 4.4 if a transposed 

765 kV horizontal line is to be modeled. The detailed design of a 
suitable L, for this t.l.m. is being carried out by Mr. William Feero 
as part of his thesis to be completed in September, 1971. Thus, the 
design of Ls is being carried out in parallel with (but slightly 
behind) the design of the remainder of the t.1.m. 

In preparation for the design of the model inductors, 
excluding consideration of the 1500 kV lines of Tables 4.4a and 4.4b, 
it was decided that the inductance values shown in Table 4.5 should 
be provided in the 4-mile sections in order to model most EHV 


lines. 
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Table 4.5 


Self Inductance Designed Into the 4-Mile Sections 


“a aD) 28 
2.24 mh 2.16 mh 2.08 mh 
2.29 ea, Yap WE) 
2.34 2.26 2.18 
2.39 2.31 DS 
2.44 2.36 2.28 
2.49 2.4] 2255 


Early in the design of the t.].m. the 1500 kV lines were 
generally ignored as beyond the modeling capabilities attainable by 
this design. With respect to the quality factor, the preceeding 
Statement is, indeed, true. However, in all other respects it is 
possible to model these lines also-with a proper choice of scale 
factor. Note in Table 4.5, under Li. the value 2.26 mh originally 
provided for modeling the typical 765 kV line. Equation (4-3) shows 
the L540 (unscaled L.) of a 345 kV line being scaled to the 2.26 mh 
inductor previously provided. The 2.26 mh value is not unique, and 
any inductance parameter of a line to be modeled could be scaled 
to any one of the inductances provided in Table 4.5. In fact, to 


properly fit the two 1500 kV lines, the L, of the horizontal 1500 


b 
KV line was scaled to 2.16 mh and the triangular to 2.41 mh. 

(The extremes of the existing L. were used in order to keep the 
resulting values of (Lo )y coo and (Li) isq9 within "reach" of the 


ranges provided in Table 4.5.) The above procedure resulted in 


an impedance scale ratio of 0.414 for the triangular line and 0.384 
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for the horizontal. The resulting power bases to be applied to the 
actual lines are 448 MVA for the triangular, and 417 MVA for the 
horizontal. Thus, the three scaling examples provided in this 
section illustrate the considerable flexibility of the "Floating 
Power Base" technique. The values of self inductance to be designed 
into the 8-mile sections are obtained by simply doubling the values 
in Table 4.5. 
4.2.2 Selection of Model Mutual Inductance Parameters 

Tables 4.4a and 4.4b indicate that typical values of Lob 
range from 0.006 mh/4 mi to 0.4 mh/4 mi, and those for Lic from 
-0.06 mh/4 mi to 0.27 mh/4 mi. Providing this wide range of values 
is neither practical nor necessary for the following reasons: 

1. To cover the range implied above would require too 

many taps on the mutual windings and too many values 

brought to the front panel of the section. As few 

as six "tap values" per component in Table 4.4 means 

that 120 lines must be brought to the front panel. 

(It is proposed to limit the number of taps on any 

component to six.) 

2. In those cases where the mutual inductance value 

is 1/20 or less of the self inductance upon which it 

is wound, it would be nearly impossible to measure 

the effect of the mutual inductance. Therefore, 

those values of mutual inductance amounting to 1/20 

or less of the self inductance in the same line will 

simply be ignored. Note that this rule leaves the 


negative mutuals ignored. 
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3. Discounting the extremes mentioned above, most 

mutual values appear to be in the range of from 

0.17 mh to about 0.27 mh - the single exception 

being the far removed 0.40 mh of the hypothetical 

type "39F" line. Mutual taps actually provided in 

the model section should reflect this “crowding” in 

the mid-band by the provision of relatively more 

taps in this band. 

4, Horizontal lines are common, and it is seen that 

for horizontal lines the mutual inductances are equal. 

The selection of taps to be provided in the model 

Section should preserve this balance if possible. 

Based upon the arguments above, the values of mutual inductance 
to be designed into the 4-mile model sections are shown in Table 


4,6 - those for the 8-mile sections are shown in Table 5.5. 


Table 4.6 
Lop mh) Lic (mn) 
eho =o 
18 18 
Al 520) 
.23 es 
BA Woe 
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Table of Mutual Inductance Values 
Designed into the 4-Mile Sections 
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4.2.3 Selection of Model Resistance Parameters 

The resistance ranges implied for a 4-mile model section are 
Shown in Tables 4.4a and 4.4b in which the figures represent the 
total resistance required of a section from input terminal to output 
terminal. Thus, the resistance indicated includes the resistance of 
the self and mutual inductance windings discussed above. The amount 
provided for operator adjustment is the difference between winding 
and connection panietance and the values indicated in the tables. 
The resistance of the self and mutual windings must obviously be less 
than the resistance values in Tables 4.4 if the 4-mile model section 
is to provide a correct modeling of line Q. The exact value of the 
resistances to be designed into the 4-mile section cannot be deter- 
mined until the resistance of the self and mutual windings is known. 
The setting of the operator adjusted resistance parameters to be 
provided is deferred to Chapter 7. 
4.2.4 Selection of Model Capacitance Parameters 

The capacitance parameters implied for a 4-mile section are 
also indicated in Tables 4.4a and 4.4b. With a wide range of 
Capacitor sizes readily available from manufacturers, a bank of 
S$1X capacitors will be provided for each capacitance parameter for 
a total of 36 capacitors in each section. The capacitors of each 
bank are to be sized so as to span the ranges indicated with "room" 
on either side. With six capacitors per capacitance parameter, 
about 4: values can be presented to cover the range of the parameter. 
Tables of the actual capacitor banks purchased for each parameter 


are presented in Chapter 6. 
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4.3 Summary Remarks 


1. The parameters of the six non-1500 kV lines indicated 
in Tables 4.1, 4.2, and 4.3 are taken to "span" the range 
of typical EHV lines (including 500 kV lines). 

2. No attempt is made to model positive sequence line 
Q's greater than 20. 

3. Six model values spanning the range of EHV line 
configuration "types" are provided for each parameter. 
4, Only three scaling techniques were illustrated, but 
any line to be modeled could be scaled to any one of the 
existing model values in an effort to make diverse lines 
fit the model. 

5. A scale factor, F, obtained for inductance and 
resistance scaling is used as 1/F for capacitance 
scaling. 

6. Values of mutual inductance less than five percent 

of the self inductance in the same line are ignored. 

7. Parameters indicated in Tables 4.4a and 4.4b for a 
4-mile model section are simply doubled to obtain the 
values for an 8-mile section. 

8. The details of the return path inductor design, Li» 
will appear in a subsequent report (by a different author) 


due in September, 1971. 
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CHAPTER 5 
INDUCTOR DESIGN 

9.1 Model Inductor Requirements 

The niodel inductance values required of a 4-mile section are 
Shown in Tables 4.5 and 4.6. The inductance values for an eight-mile 
Section are simply double those of Tables 4.5 and 4.6. In addition 
to simply providing the inductance, the inductors in the model sections 
lust have Q. considerably greater than the Q characteristic for the 
"best" line to be modeled. The Q characteristic of the 765 kV 
norizontal line shown in Figure 4.1 is taken as the minimum standard 
to be met by the model sections. The Q characteristics of the 
individual model inductances inust exceed the requirements of the over- 
all section. while the line and section require a Q of twenty at 
00 Hz, the inductors should have a Q of perhaps 40 at 60HzZ - anda 
Similar excess over the entire frequency range. The characteristic 
shown in Figure 4.1 is the “average’ Q of the three phases so that 
it is not the exact characteristic to be tracked by any of the 
three self inductors of the model even if this tracking were possible. 
Tne solution of the "Q Tracking Problem" in this design is to design or 
Select an inductor whose Q characteristic has the general shape shown 
in Figure 4.1 but whose Q value is approximately fifteen to twenty- 
five points higher along the range. Any excess Q of the model can 
be handled by simply inserting an additional small resistance in 
series with the model inductors. 

In addition to the Q characteristic per se, the model 


inductors inust produce the required Q at the very low values of base 
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current selected for the niodel, and must have an inductance that remains 
essentially constant up to at least 25 p.u. fault current. The 
inductance of all the phase and mutual inductors iust also renain 
constant with frequency - only the return path inductor, L 1s to 

vary with frequency to provide the zero sequence inductance characteristic. 
9.2 Inductor Core Selection and Core Properties 

Tne inductor characteristics described in section 5.1 are 
deteriiined, to a very large extent, by the core material of the inductor. 
Tne theory leading to tne final relations describing inductor behavior 
is rather involved and will not be included in this report - only the 
More significant relations will ve presented to indicate the overal] 
procedure. Readers interested in a more complete background leading to 
tne relations presented herein should consult references 13-19. 

A number of core materials and configurations were investigated 
to meet the specifications described above, but most were readily found 
to be unsuitable and rejected early in the design. A potted ferrite 
core, for example, meeting all other requirements, was rejected as too 
expensive (Ferroxcube series 6656). Two alternatives were considered 
in some detail - oriented silicon steel "C" cores and molybdenum- 
nickel-iron powdered toroids. The oriented silicon was rejected when 
the model base current was finally set at the relatively low value of 
38.6 tia. At this low current the oriented silicon could not produce 
the desired Q characteristic at a sufficiently low cost or weight. 
Basically, the Q of oriented Fe-Si inductors increases with increasing 
current (up to saturation), whereas with molybdenum-nickel, Q 
decreases with increasing current. Tnus, the Mo-Ni cores nave their 


highest Q's at low currents and the Fe-Si cores have their lowest Q's 
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at low currents. The basic relation governing the Q of an inductor can 


be written as (for the frequencies of interest here) 


2 
= Wh. ; (5-1) 
Rael 1: Racl 
Where 
L = Inductance in henries 
w = 2m x frequency 


Rac = dc resistance of the winding (2) 


R 


_ equivalent core loss resistance (&) 


With these definitions it is seen that R, oI corresponds to both the 
hysteresis and eddy current losses (See reference 17 p.p. 58ff.) and 


thus represents the total core loss. Rae? the core loss coefficient, 


is itself a function of current. For powdered cores the manufacturer! © 


provides data on se explicitly, whereas for Fe-Si the entire term 


RI is provided. If the data provided by the manufacturer 1s 


converted to the terms of equation (5-1), it is found that for Fe-Si 
the core loss, W, is given by (approximately) 


We ki (wll? (5-2) 


] 
whereas for the powdered Mo-Ni the loss is given by (approximately) 


W=R 1° 3 ko (w) 1% 4 : 


ac (5-3) 


In equations (5-2) and (5-3), ky (w) and ko(w) are core 10Ss constants 
and are a function of frequency. If equations (5-2) and (5-3) are 


Substituted in equation (5-1) there results, respectively, 
2 
= -. obloy 2 (5-4) 


Q ~ e 
oa D ey, 
del + ky (w)I 


and 
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wh? 


2 TT (9) 
Rac! + ko (w)I 


Soetlia 


Equations (5-4) and (5-5) apply only in the region below saturation 
and suggest the reason wily the Q of the Fe-Si core increases with 
current, whereas in the Mo-Ni core Q decreases. The Q calculations 
were conducted for a number of oriented Fe-Si cores (Arnold-AH-393, 
AH-34, AL-59, AH-95, AH-86, AH-31), maximized with respect to air gap, 
with the result that cores up to eight pounds each, costing up to 
$39.00 each, could not produce the required Q characteristic at the 
base current levels required. 

Several Mo-Ni cores could meet the characteristic, but the 
Smallest (in terns of those offered by the Arnold Engineering Co. or 
ilagnetics Corp.) and least expensive unit which could be made to 
produce the desired ciiaracteristic with a Q of 40 at 60 Hz was Arnold 
part no. A-218259-2. The Q characteristic of this small (2.25 inches 
Q.D.) toroid when full wound is shown in Figure 5.1. The curve above 
100 Hz is the manufacturers specification which was verified for a 
typical full wound core using a GR 1650A inductance bridge down to 
300 Hz. The curve below 1U0 Hz is the extrapolation of the author. 
The Q at 59 Hz was checked using the circuit shown in Figure 5.2. 
R. was determined by measuring the dissipation factor, D, of the 


870 uf capacitor bank and using the relation 


rp =2 (5-8) 


Coen 
D was measured at 60 Hz using a GR 1611-A capacitor bridge and found 
to be .0045. Five copper lug connections are included in the tank 


circuit along with fhirteen feet of No. 10 AWG solid copper magnet 
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wire for a total measured resistance of 15 mQ using a GR 1650-A 
resistance bridge. Ine voltage waveform across the inductors was 
observed witn a Tektronix 549 storage oscilloscope. The 8.4 mh 
inductor Shown in the figure is four 2.1 mh inductors wound on the 
t.].m. cores. The tank circuit was "pumped up" using the oscillator 
at 59 Hz after which the switch was opened and the decay waveform 
observed on the Storage oscilloscope. Twelve separate measurements 
were miade covering the range of initial voltages from 0.02 v peak 
to 0.4 v peak - corresponding to peak initial currents in the range 
from 7 ma to 129 ma. In each case the time was measured for the 
voltage across the tank to fall tol of its initial value. The 
times measured ranged from a low of 141 m sec at 0.4 v peak to 

157 m sec below 0.1 v volt peak. The total resistance in the tank 
was detennined frai 


_ 2b _ 16.8 mh 


Rt = time ~ “tine — ) 
Where | 
Bye herd Reet Rae Roe (5-8) 
Finally, Q of the inductors was determined froin 
L 
Q = We (5-9) 
Rac i Rac 


The twelve Q values at 59 Hz ranged from a low of 34.6 at .4 v peak 
(initial I = 129 ma), to a high of 40 for initial voltage values 
below 0.1 v peak (32 ma). A similar “pumping" technique was used 

to measure the inter-winding capacitance of a single t.].m. inductor 
of 2.1] mh. A series circuit was formed of an H.P. 4204A variable 


frequency oscillator, a 4.7 meg 2 resistor, and a Single 2.1] mh 
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t.].m. inductor. The voltage across the coil was observed while the 
frequency waS increased. The voltage peaked sharply at 300.0 KHz 
implying an inter-winding capacitance of 134 p.f. 

The maximum inductor Q shown in Figure 5.1 occurs when the 
losses in the core of the inductor equal the losses in the resistance 
of the windings. At low frequencies the winding losses predominate; 
at maximum Q the two losses are equal, and for frequencies above 


Q 


Q curve is determined by winding losses which remain essential ly 
17 


a the core losses predominate. The initial slope of the inductor 


constant in this frequency range. In other words - at low frequencies, 
the Q is determined by winding resistance alone. In fact, a t.].m. 
inductor of 2.1] mh at 60 Hz has wl = 800 m2 and a winding resistance 
of 20-m2 indicating a Q of 40. To bring the Q of the t.].m. section 
down to 20, for example, 20 m2 of resistance is to be inserted into 
the inductor winding - sufficient to bring the overall Q of the 
section at 60 Hz down to about 20. Other, larger values of resistance 
would have to be added to produce Q's less than 20 for modeling lower 
Q lines. The second curve in Figure 5.1] represents a typical Q 
characteristic of a section. Note that the section characteristic 

is lower, and the maximum is shifted to the right - reflecting the 
fact that core losses must reach a higher value before they can 

equal the winding losses thus producing the Oak condition. Note 
also that having shown a case in which the winding resistance was 
doubled, the initial slope of the section Q curve is about one-half 
of that of the inductor alone. The third curve in Figure 5.1] is that 


of the typical 765 kV line transferred from Figure 4.1]; this curve 
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represents the highest Q which must be achieved by a section. It is 
seen that setting the Q of the section at 60 Hz to twenty results 
in the section Q falling below that of the line above 1700 Hz. Thus, 
if Q is a critical factor above 1700 Hz in a particular test on a 
typical 7/65 kV line, it will be necessary to accept higher Q at sixty 
Hz by setting less than 20 m2 in series with the model inductors. 
9.3 Design of Ilodel Self Inductors 

The inductance of a toroid is given by 


ane uA x 107° 
L = ———~——_—_—— henries, (5-10) 


D 
where 
N = number of turns 
= relative permeability 


2 
core cross sectional area (cm) 


> F 
Ul 


D = mean core diameter (cm) . 
For the Arnold core A-218259-2 used in this design, uM, = 205, so that 
L = 259 mh per 1000 turns on any full wound core. A full wound core 
is defined as one in which the winding occupies three-fourths of the 
available window area with a packing factor of 0.7. For any other 


full wound coil on the A-218259-2, 


: J. : 
N = 1000 /seq ; (5-11) 


where "L" is the desired inductance in mh. The inductances shown in 
Table 4.5 require the turns shown in Table 5.1. The listing of turns 
required in the 8-mile sections is shown in Table 5.2. To achieve 


the Q characteristic shown in Figure 5.1, it is necessary to have the 
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Table 5.1 


Turns Required in 4-Mile t.1].m. Inductors 


"a wa Nc 
93 9] 90 
94 a2 ot 
95 93 32 
96 94 g3 
97 95 94 
98 96 sls 
Table 5.2 


Turns Required in 8-Mile t.1].m. Inductors 


Wa "> 2 
13] 129 ee 
SIS: 130 129 
134 13] 130 
136 133 Vl 
137 134 133 
139 136 134 


core full wound. In the case of the 4-mile inductors this works out 
to be No. 10 AWG, and for the 8-mile sections - No. 11 AWG. For the 
4-mile sections, a nominal twenty feet of No. 10 is required at 1.0 
m2/ft for 20 m2 per inductor. For the 8 mile sections, a nominal 
twenty-eight feet of No. 11 is required at 1.26 m2/ft for 35 m 


per 8-mile inductor. In Tables 5.1 and 5.2, the first five numbers 





8/7 


in each column are the turn numbers at which a tap must be provided. 


Referring to equation (5-10) it is seen that the inductance 


of a toroid is independent of current if HW is independent of 


current. A graph of H, VS. a-c flux density taken from reference 18 


is reproduced in Figure 5.3. 


percent of nominal up to 4000 gauss. 


b Heds = NI 
TDH = NI 
= NI/nD 


Substitute (5-15) into (5-12): 


NBA 
I 


=| CO 


I 


In equation (5-18) 
I 
N 
D 
B 


u 


Ne uA 


TD 


complex amplitude of current (amperes ) 
no. of turns in winding 

mean diameter of core in meters = .0463 m 
peak flux density (wb /m) 


normal permeability = up 


For the t.l.m. inductors at 4000 gauss: 


B 


u 
D 


N 


r 0° 
4 wb /m¢ 
205(4nx107/)= = 258x107" 
0.0463 m 


3703 


From basic inductor theory: 


It is seen that Mi remains within five 


(5-12) 


(5-13) 
(5-14) 
(5-15) 


(5-16) 


(5-17) 


(5-18) 
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m™DB 
T33u amperes 


So that I peak 


= 1.70 amperes, 

and the inductance of the t.].m. inductors remains within five percent 
of nominal up to 1.70 amperes r.m.s., or 31 p.u. 
5.4 Design of Model Mutual Inductors 

Figure 2.2 shows the mutual inductors Ly o(L.) and L53(b..) 
as transformers between the phases. Since the mutual inductances are 
small, and since there is sufficient space remaining in the windows 
of the self inductors, L eb is wound over Lo and Lie 1s wound over Lh 
The mutual inductances to be provided in the 4-mile sections are 
shown in Table 4.6; those for the 8-mile sections are simply double 
the values in Table 4.6. From basic mutual inductor theory (See e.g. 


reference 15 p.p. 178 ff) it is known that, for example, 





Lob = k /(Lon). L. (5-19) 
ee ae 7 
be aati i (5-20) 
and from equation (5-11) for the model core, 
(N_.). = 1000 Mab’ (5-21) 
ab’s 259 : 


Where (all in mh): 
Lob = mutual inductance between phases "a" and "b" 
(en es self inductance of L.,. 
k = coefficient of coupling between L. and Lb 
L" = measured mutual inductance when L. and Lb 
are connected to reinforce. 
L" = measured mutual inductance when L, and L., 


are connected to oppose. 
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(Nand. = No. of turns to obtain (L or E apie 


ab/s* ( 
Several tests were performed to determine k by winding six, seven, and 


eight turns of (Li) over an L, inductor of 2.49 mh. Knowing the 


S 
number of turns, (Lins was obtained from (5-21), and Lob was obtained 
as indicated by (5-20). ae and Lob were then substituted into 
(5-19) to determine k. Under these conditions, k for the t.1.m. 
inductors was found to be 0.986. Solving (5-19) for (Lind. and 
Substituting into (5-21), 


1 


Oe = 63 Lob tL : (5-21) 
With L. taken as 2.49 mh, 
ieee = 39.8 Lob? (5-22) 


and the number of turns to achieve the mutual inductance values, Lb? 
of Table 4.6 are 6, 7, 8, 9, 10, and 13, so that taps must be provided 
at 6, 7, 8, 9, and 10 turns. With the taps set, the self inductance 
corresponding to each tap can be found using equation (5-21) and then 
the values of Lib with other values of L, can be found using (5-19); 
the results are shown in Table 5.3. A similar procedure applied to 


L.. yields, 


bec 


(th. Gldois I 


be be 

and still the taps are at 6, 7, 8, 9, and 10 turns, with the 13th 
turn providing Lic = 0.32 mh. A list of Lic values for a 4-mile 
section is given in Table 5.4. In order to wind Lb and Lic over 

L. and Li. respectively, it is necessary to use four feet of No. 12 
AWG magnet wire at 1.59 mQ/ft for an additional resistance of 6.4 m2 
in series with the phase "b" and "c" lines. A list of mutual 
inductance values to be built into the 8-mile sections is shown in 


f 
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Table 5.5. 
Table 5.5 
Lop wma ) Lic kmh) 
0.30 0.30 
0.36 0.36 
0.40 0.40 
0.46 0.46 
0750 0-50 
0.66 0.56 


Table of Mutual Inductance Values 


Designed into the 8-Mile Sections 


The taps to yield (Laide (using La 4.98) are 8, 10, 11, 13,1425 and 
18; taps for (Li de (using Peas 4-82)ware Siel0, Wik 13.84 ewandel7: 
The corresponding values of (Ln). for the 8-mile sections are (in mh): 
Wmo1o6, 0.0259, 0.0312, 0.0435, 0.0505, and 0.0834. For (Li de: 
0.0166, 0.0259, 0.0312, 0.0435, 0.0505, and 0.0745. These tabulations, 
along with those of Tables 5.3 and 5.4 show that the self inductance 
of a mutual winding is negligible when compared to the self inductance 
of the line. In addition, Tables 5.3 and 5.4 indicate that the change 
in ah or Lie resulting from a change in L. or Ly over their entire 
range is less than five percent and therefore negligible also. In 
order to wind Lab and Lie over L. and L. in the 8-mile sections, it 

is necessary to use approximately six feet of No. 12 AWG magnet wire 
at 1.59 m2/ft for an additional resistance of 9.5 m2 in series with 


the phase "b" and "c" lines. 
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CHAPTER 6 
CAPACITOR SELECTION 

6.1 Factors Determining Capacitor Selection 

As discussed in detail in Chapter 4, eight diverse transmission 
line types were analyzed in the EPSEL Line Parameters Program, and 
these eight types were taken to cover the range of Typical EHV lines. 
The model capacitance values implied by these lines for a 4-mile model 
section are shown in Tables 4.4a and 4.4b. A complete listing of the 
capacitance parameters to be covered by both the 4 and 8-mile model 
sections is shown in Table 6.1. Having scaled the eight representative 
lines to result in a minimum inductance range, the range of each 
Capacitance is seen to be somewhat larger. It is proposed to cover 
the range of each parameter of Table 6.1 with a bank of six capacitors 
wired and brought to the operating panel so that any combination of the 
Six capacitors in each bank can be arranged in parallel. Thus, with 
one or two capacitors to set the lowest value in a range, approximately 
five or four capacitors are available to cover the remainder of the 
range. Thus there are a maximum of 120 and minimum of 24 capacitance 
values available to cover the range of each capacitance parameter 
indicated in Table 6.1. With this "dense" provision of model values 
it is possible to set any value in the span within a very few percent 
(v3%) if the tolerance of the capacitors is less than a few percent. 
It is proposed to provide a two to three percent tolerance on any 
capacitance selected by the operator so that the over-all capacitance 
provision of the model is within about 5% of the value required to 


model any given line. With expected transient overvoltage maxima 
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of about four p.u., all capacitors must be capable of withstanding short 
duration peaks up to 500 or 600 volts. In addition to tolerance and 
voltage requirenents tie capacitors must have extremely low losses. 

An indication of tne dissipation factor required of the model 
capacitors can be obtained by a simple analysis of the "dissipation 
factor" of a typical line. Considering the positive sequence parameters 


of a type 41A line using the values given in Table 4.1, 


B 
cee 2s 0181 uf/mi. 
Un = ,U/2o0t/4. ii. 
Ros. = .216 2/4 mi. 
Using the "series model" for a lossy capacitor: 
Ups /4 mi = WRp g Cp ig. (6-1) 
Up S = 0.0000059/4 mi. Q 60 Hz (6-2) 
Up 5 F 0.0000985/4 mi. @ 1900 Hz (6-3) 


That the values given by equations (6-2) and (6-3) are the values 
required of a 4-mile model section can be seen by a corresponding 


calculation for the model. Using Table 4.4a for the type 41A model: 


oe ae LO + C 
me! T _ “an bn cn ; 
pasa. <0) cei neem O Ra ace 
CDs = .230 uf/4 mi. 
Rk +R +R 
T _osoa b C _ . 
Rp = Re = —— oa = .068 2/4mi. 
Up SS 0.0000059/4 mi. @ 60 Hz (6-4) 
8 = 0.0000985/4 mi. Q@ 1000 Hz. (6-5) 


Boe 
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at low currents. The basic relation governing the Q of an inductor can 


be written as (for the frequencies of interest here) 


2 
Q ee eee (5-1) 
Rae! Ract 
Where 
L = Inductance in henries 
w = 2m x frequency 
Rac = dc resistance of the winding (&) 
Neve = equivalent core loss resistance (Q) 


With these definitions it is seen that R, oI corresponds to both the 
hysteresis and eddy current losses (See reference 17 p.p. 58ff.) and 


thus represents the total core loss. R the core loss coefficient, 


ac’ 


is itself a function of current. For powdered cores the manufacturer! © 


provides data on ee explicitly, whereas for Fe-Si the entire term 


R If 


- is provided. If the data provided by the manufacturer is 


converted to the terms of equation (5-1), it is found that for Fe-Si 


the core loss, W, is given by (approximately) 


We ky(w)d'*7, (5-2) 
whereas for the powdered Mo-Ni the loss is given by (approximately) 
: 2 2.4 : 
aL ee ea AC )) (5-3) 


In equations (5-2) and (5-3), ky (w) and k5(w) are core loss constants 
and are a function of frequency. If equations (5-2) and (5-3) are 


Substituted in equation (5-1) there results, respectively, 


1. coe wh ® (5-4) 
Fe-S] Ryol’ ri k (wt! 


and 
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occurs primarily at low frequencies in connection with switching opera- 
tions which connect the terminals of a charged line to large inductive 
reactances. Natural frequencies of the order of 40 Hz are known to occur 
under these conditions with over-all Q approaching 300. 

Most capacitor manufacturers specify or test dissipation factor at 
1000 Hz, so that at 1000 Hz it is possible to measure capacitance and 
compute resistance using either the parallel or series model of a lossy 
capacitor: 


7 | 
Dowie Poe . (6-6) 


Pe 


The values of R. or R_ so obtained, however, are true only for a limited 


p 
range of frequencies and cannot be used at 40 Hz. As frequency decreases 
toward zero the value of Re (resistance in the parallel model) approaches 
the insulation resistance of the dielectric - which for polystyrene is 


20 The 


of the order of ten million megohms for a 0.15 uf capacitor. 
dissipation factor, D, of a typical .15 uf polystyrene capacitor at 

1000 Hz is given as 0.0001.°? This value of D implies R, = 10.6 megohms 
at 1000 Hz. Few materials other than Teflon or a modified polystyrene 


can equal this low value of D or high value of Ro ac l000 iiz. min 


equation (6-6) is solved for Re in terms of w, Re, Cs and C. 

~ 
(C. = Ca) and Re set equal to 1.0 m2, C. = 0.15 ut, and | = 40°Hz the 
value of R. obtained is 4.4 million megohms. Thus in order to obtain 


p 
ie as low as .00] 2, it is necessary to have Ro as high as 4.4 million 


megohms at 40 Hz. Little data exists on any dielectric material below 


100 Hz - standard reference works<- 


supply no data below 100 Hz on any 
material and give no information vs. frequency on polystyrene. In an 


effort to obtain an estimate of Re for the 0.15 uf polystyrene capacitor 





he, 


at 40 Hz, it was assumed that Ro remains constant at the 10.6 megohm 
value down to 100 Hz and then rises linearly to 10 million megohms at 


d.c. This procedure gives R. as 6 million megohms at 40 Hz. (This is 


p 
a conservative estimate in the opinion of the author.) No material 
less costly than polystyrene can approach within an order of magnitude 
of this value. Polystyrene capacitors are therefore selected to form 
the “base” capacitance of each capacitor bank in the model. 
6.2 Capacitors Selected for Model Sections 

Since the polystyrene capacitors conservativety (but not by 
orders of magnitude) meet all line requirements, a compromise in favor 
of lower cost is made in selecting the trimming capacitors in each 
bank. Thus, while polystyrene capacitors will be made to provide 
80 percent (approximately) of any capacitor value selectable by the 
Operator, the remaining 20 percent (approximately) is to be made up of 
lower cost mylar or silvered mica capacitors. The latter capacitors 
have losses ten to fifty times higher than the polystyrene but they 
constitute only ten to twenty percent of a selected capacitor value 
and so should not appreciably lower the high quality provided by the 
polystyrene base capacitors. The severe low-loss constraint applies 
mainly to the 765 kV line, so that the polystyrene capacitors in each 
bank are selected so that for the 765 kV case they will constitute no 
less than 90% of the required capacitance of a 4-mile section. The 
polystyrene capacitors selected for the t.].m. are Arco Manufacturing 
Co., type PJ, 400 VDC, + 2%, tubular capacitors. The 400 V units 
were judged adequate since each capacitor receives a minimum of two 
sixteen hour burn-in periods with 300 percent of rated voltage applied 


20 


aleeco C, The mylar capacitors are Cornell Dubilier Co., type WHF6, 
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600 VDC, + 10 percent, tubular capacitors. Jwo small capacitance values 
below the range of standard mylar sizes (330 pf and 680 pf) are Cornell 

Dubilier Co., type CDI9FD, silvered mica, 500 VDC, + 5%. The capacitor 

types and sizes selected to span each capacitance parameter (the 


capacitor banks) are shown in Table 6.2. 
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CHAPTER 7 
RESISTANCE AND LOSS CONSIDERATIONS 

7.1 Over-all Q Considerations 

As discussed in Chapters 4 and 5, the Q requirements of this 
t.1.m. are dictated by the need to model the 765 kV line. This con- 
Straint is particularly severe in the 4-mile model sections since 
connection resistance is expected to take up any margin provided by 
designing the inductors with double the required Q at 60 Hz. In 
Chapter 4 it is indicated that the positive sequence Q of a typical 
765 kV line is twenty. The resistance required in each phase to achieve 
this result is shown in Table 4.4a - 37 m2 in phases “a" and "c" and 
43 m2 in phase “b'. As discussed in Chapter 5, 20 m2 is "provided" in 
each phase by the resistance of the wire to wind the self inductors 
Lo Lie and Lee An additional 6 m2 is provided in phases “b" and “c" 
due to the fact that these phases contain the resistance of the windings 
for L.. and Lic respectively. (Actually, Lab and Lac for the typical 
765 kV line are provided by the 7-turn-tap so that the resistance 
added by the Lb and Lie windings will be approximately 4 mQ in the 
4-mile sections. The operating panel resistance to be provided at this 
point in the accounting is 17 m2 in phase “a", 19 m2 in phase "b", and 
13 m2 in phase "c". But connection resistance has not yet been 
accounted for in these figures. 

In a typical model set up (See Chapter 8), six to eight banana 
jack connections will be made in series in each section from input to 
output, so that the resistance of these must be accounted for. The 


jacks selected for this design are E.F. Johnson Co., standard banana 
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jack, cadmium plated insert, catalog no. 108-0900-001. The corresponding 
plug is E.F. Johnson Co., standard plug, catalog no. 108-0300-001. Six 
of these jacks and plugs were connected in series as shown in Figure 7.1. 
Sixteen inches of solid, no. 14 AWG, copper bus bar wire was used, and 
all connections were heavily soldered. The total resistance was measured 
at 12 m2 with the GR 1650A resistance bridge. The resistance of the 

no. 14 wire is 3.3 m2 so that the resistance of each connection averages 
to about 1.5 m2. All plugs and jacks were new, clean, tight fitting, 

and all connections were made and remade several times - always with the 
Same result. The resistance accounting to this point is: 


R, = 8 m2, R 


, a 5 mo, R. = 1] m2. In other words, at this point the 


resistance values required at the front panel to bring the section 
positive sequence Q down to twenty at 60 Hz are the values cited just 
above. Each 4-mile section; exclusive of external hookup wire, but 


including operator panel connector resistance is seen to have built 


into it, 
R (intrinsic), = 29 m0 (7-1) 
nV) 
R.(intrinsic), = 38 m (7-2) 
a e e e V 
R (intrinsic), = 38 mQ (7-3) 


To provide an exact positive sequence Q of twenty at 60 Hz, it is 
necessary to provide for operator selection; RS = 8 mQ, R, = 5 m, 

and R = 1 m2. As pointed out in Chapter 5 in connection with 

Figure 5.1 however; setting the Q to twenty at 60 Hz for a 765 kV line 
model will result in the model Q falling below line Q above about 


1700 Hz. Thus, it may be desirable under certain circumstances to set 
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Test Circuit for Measuring Resistance 
of Banana Jacks and Plugs 


Figure 7.1] 
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Q at 60 Hz slightly above twenty to ensure that the model Q remains 
above line Q up to 3060 Hz. 
An accounting procedure identical to that above applied to 


the 8-mile model sections gives the results, 


R, (intrinsic), = 44 mo (7-4) 
e e e AU 

R (intrinsic), = 57 mQ (7-5) 

R (intrinsic), = 57 mo (7-6) 


The values of Ra? Ry? and R required in the 8-mile sections to model a 


765 kV line are (doubling the appropriate values in Table 4.4a): 


(Ree = 74 mo Gi) 
ae = BE me (7-8) 
(ere = 74 mo (7-9) 


It is seen that the 8-mile sections have considerably more margin than 
the 4-mile sections and, in fact, using the intrinsic values of 
equations (7-4)-(7-6), the positive sequence Q of the 8-mile sections 
can be set as high as 31 at 60 Hz. 
7.2 Resistance Selection 

Referring to Table 4.2 (excluding the 1500 kV lines) the values 
of positive sequence Q to be provided are seen to range from about 
eight to twenty. To achieve the Q values lower than twenty, resistance 
values must be provided at the operating panel for insertion with the 
intrinsic resistance of the section. For the 4-mile sections, the 
lowest resistance value to be made available is clearly zero since 
this is the required value to ensure positive-sequence Q's just over 
twenty. For the 8-mile sections a similar zero base will be provided 


Since it would simply not make sense to deny the operator the selection 
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of positive sequences Q's to thirty should the occasion arise. A list- 


ing of Q values to be provided for each phase is shown in Table 7-1. 


Table 7.1 
4-Mile Section 8-Hile Section 
Qs Q Q Qa Qb Q 


30 22 21 40 30 28 
Go 19 18 34 a 24 
20 16 15 28 20 20 
15 13 12 22 16 16 
12 1] 10 16 12 1] 
10 g 8 1] 9 8 


Table of Q Values Provided 
in the T.L.M. Sections 


The operating panel resistance values to achieve the Q's shown in 


Table 7.1 are shown in Table 7.2. 


Table 7.2 
4-|lile Sections 8-llile Sections 
Rs Rb Ro RS Rb Ro 
0 0 0 0 0 0 
6 7 8 8 1] 12 


15 15 17 19 28 25 
30 2/ 30 36 49 46 
44 40 43 66 85 92 
Sy) 67 64 116 132 147 


Table of Operating Panel Resistance Values 
Corresponding to Q's of Table 7.1. (R in mQ) 
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The resistance values indicated in Table 7.2 were obtained using 


w(L. ) 
pe eee anerinaies (7-10) 
1 Q; 1 
(Lida = 2.34 mh, (Li), = 2,26 mh, (Li), = 2,18 mh 
(L.)e = 4,68 mh, (Lids = 4,52 mh, (Lode = 4,36 mh 


1 = a,b,c for each section length 


Q, is the corresponding Q value from Table 7.1. 


It is seen in practice, rather than setting positive sequence ( per se, 
the Q of each phase is set individually using Tables 4.4, 7.1, and 7.2. 
The resistance indicated in Table 7.2 is achieved by connecting 
resistances in series so that each indicated value is the sum of the 
resistances required to form each of the preceeding resistors. The 
actual resistors required to obtain the operating panel values indicated 


in Table 7.2 are shown in Table 7.3. 


Table 7.3 
4-‘lile Sections 8-Mile Sections 
K S Re S RO S K S Ry S Re S 
6 7 8 8 1] 2 
9 8 e) 1] 17 ig 


15 iz 13 17 2) 2] 
14 13 13 30 36 46 
15 2/ 21 90 47 59 


Table of Resistors in the T.L.M. Sections 


The resistors indicated in Table 7.3 are to be "manufactured" by 


cutting appropriate lengths of nickel or copper-nickel alloy ("cupron") 
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resistance wire. <A very large selection of resistance wire sizes and 
types are available from the Driver Harris Co. of Newark, fi.J., a type 
particularly suited to the smaller resistors above is the company's 
0.035" diameter nickel wire of 0.045 ohms per foot. For the larger 
resistors no. 20 AWG "Cupron" (copper-nickel) wire with a nominal 
resistance of 0.28/ ohms per foot can be used. The latter wire can 
safely carry a direct current of three amperes indefinitely but is 

hot to touch at this current level. Three rms amperes is far in 
excess of that required to simulate the most severe faults, so that 


no. 20 AWG resistance wire should be adequate for this t.l.m. 
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CHAPTER 8 
MODEL PHYSICAL CONSTRUCTION 
8.1 Electrical Construction of a Typical Section 
As discussed in Chapter 4, this t.].m. consists of thirty 4-mile 
sections and thirty-five eight mile sections. A simplified schematic 
diagram of the wiring of a typical section is shown in Figure 8.1 and 
is the same for both the 4-mile and 8-mile sections. The component 
values are not indicated since they are different for the two section 
types, but the required values can be taken from the tables of Chapters 
5, 6, and 7. In wiring each section, all hookup wire in a direct path 
from input to output is No. 10-stranded, tinned, copper wire. Capacitor 
hookup is effected with No. 16 stranded copper wire. All connections 
are soldered. 
8.2 T.L.M. Physical Construction 
Each t.].m. section is constructed on the chassis shown in 
Figure 8.2. The sixty-five chassis are made of plexiglass of the sizes 
indicated - the operator's panel is of black plexiglass. The layout 
of the operator's panel is shown in Figure 8.3. The jacks are color 
coded by "phase" and the circuit diagram is painted by a silkscreen 
Stenciling process. Each of the sixty five chassis is assembled in a 
larger console as shown in Figure 8.4. The chassis slide into the 
console on rails and are secured with screws through the operating 
panel. The chassis are inserted into the console from both "front" 
and "rear", so that twenty-four sections are accommodated per console. 
Three consoles such as that shown in Figure 8.4 are required with a 


fourth "half-console' for the output section and terminal equipment. 
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CHAPTER 9 
SUMMARY AND RECOMMENDATIONS 
9.1 Summary 

This report presents in detail the considerations in the design 
of a transmission line model for use in the M.I.T. Transmission System 
Simulator. The design includes the modeling of the inductive effects of 
non-transposition by incorporating two linear mutual inductors between 
phases. The inter-phase resistance effects are approximated by assuming 
that the three inter-phase resistances are equal. A circuit diagram of 
one three-phase L-section is shown in Figures 2.2 and 8.1. Note that 
the mutual inductors are physically wound on the respective self inductor 
core and that the return circuit is provisional pending a complete 
design due in September, 1971. 

In order to deduce the number of sections required for the model 
and establish a quantitative measure for predicting the performance of 
model types with various terminations, the pole error analysis of Chapter 
3 was developed. This analysis is based upon a comparison of the poles 
of the driving point impedance of actual lines with those of proposed 
models. Solutions for the poles of 7 and L-section models were obtained 
by computer for models up to forty sections with open and short circuited 
terminations. The pole error curves of Chapter 3 show in quantitative 
terms the convergence of the poles of various models toward those of 
continuous lines with similar terminations. A technique is outlined 
(See Figure 3.11) in which the pole error curves are used to deduce 
the number of sections required to meet any arbitrary accuracy spec- 


ification. The specification set for the M.1.T. model is that the 
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error in the pole frequency of the model shall be within five percent 
of the frequency of the corresponding actual pole over the range from 
sixty to 3060 Hz for line lengths of sixty to four hundred miles for 
both open and short circuited terminations. This criteria led to the 
decision to construct thirty four-mile model sections and thirty-five 
eight-mile sections ("L" configuration). 

In Chapter 4 parameter data for eight representative EHV line 
configuration types are presented in order to deduce the parameter values 
required in the four and eight-mile model sections. Several scaling 
techniques are suggested for scaling the parameters of the diverse line 
types to the parameters provided in the model. Careful scaling suggests 
that only a narrow range of component values need be provided for each 
parameter in the model. In the M.I.T. design, six values are provided 
for each of the sixteen parameters to be modeled. 

The quality factor (Q) of typical EHV lines are seen (Figure 4.1) 
to vary markedly with frequency in the range from sixty to 3060 Hz. 
Positive sequence Q increases monotonically to over 175 for the highest 
Q lines studied, and zero-sequence Q maximizes near 1000 Hz and 
generally remains constant up to 3060 Hz. The requirement to model 
these Q characteristics is a major factor in the design of all model 
inductors and in the selection of all model capacitors. The positive 
sequence Q characteristics is achieved by the use of molybdenum- 
nickel toroids for all phase inductors. In fact, the model inductors 
have an intrinsic Q characteristic approximately double the value 
required to model the highest Q EHV lines in operation. Preliminary 


plans for modeling the return path or zero-Sequence inductance and 
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Q suggest a separate modeling of the ground wire path and earth return 
path - the final design is to be completed by September, 1971. 

A careful account of wire and connector resistance shows that 
for the section lengths chosen, every milliohm significantly reduces 
the intrinsic Q provided by the model inductors. To avoid excessive 
losses and provide for operator adjustment of Q, all connections are 
soldered, and all wire on a direct path from input to output is 
no. 10 AWG hook up wire. To avoid further losses, and to model the 
condition of a charged line isolated with a high Q reactor at its 
terminals, extremely low loss polystyrene capacitors are incorporated 
into each t.1.m. capacitor bank. Approximately eighty-five percent 
of any mode] capacitance parameter selectable by the operator is low 
loss polystyrene in parallel] with approximately fifteen percent 
additional "trimming" capacity provided by mylar units. Overal] 
Capacitance tolerance in the t.].m. is within five percent of any value 
selected. Six capacitors are provided for each capacitance parameter 
and all capacitors are available for operator connection into any 
desired configuration. 

The final result for the M.I.T. t.1].m. is a total of sixty-five 
L-sections for modeling line lengths to 400 miles. Average Q in the 
four-mile model sections is a maximum of 24 at 60 Hz; the maximum 
average in the eight mile sections is 32. All capacitance and 
inductance selection within the range of a parameter is within five 
percent of any desired value. A similar five percent tolerance is 


planned for the return path inductance and Q. Q values less than the 
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maxima given above are obtained by the insertion of appropriate resistors 
provided in each section. 
9.2 Recommendations 

The design of this t.}].m. has progressed in parallel with its 
actual construction and component procurement. This was attempted in 
view of a master's thesis time constraint. The basic scale factor 
between the model and lines was selected in order to make the t.1.m. 
compatible with existing laboratory equipment. In fact, the t.1.m. 
ultimately will represent an investment several times the cost of the 
TSS transformer which dictated the scale ratio. This scale ratio 
(the 0.355 for the 765 kV line) was selected and found acceptable 
primarily on the basis that it resulted in a feasible, low cost inductor 
core without scaling the capacitors out of the range of economical 
values provided “off the shelf" by manufacturers. Not until after the 
cores had been purchased was the significance of capacitor dissipation 
factor fully appreciated. Although the required capacitor values are 
well within the range of normal capacitors, they are among the highest 
and most costly polystyrene units. The result was the requirement to 
purchase capacitors nearly three times the cost of inductor cores. 
Polystyrene capacitors are required regardless of what reasonable scale 
ratio is selected, but it is suggested that perhaps a better cost trade 
off could have been achieved had the TSS transformers been redesigned to 
permit a scale ratio closer to unity. This would require larger, more 
expensive, cores for the t.].m. inductors but would significantly reduce 
the size of capacitors required. The larger polystyrene units are very 


expensive, so that even a small reduction in values required could result 
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in considerable savings. This savings could be directed toward purchasing 
a better (larger) inductor core permitting even more margin in Q by making 
the connector resistance a smaller fraction of the total resistance of the 
section. It is recommended that, should a similar project be undertaken 
in the future, the entire design should be completed in every detail, a 
prototype section constructed, and a cost analysis completed before a 
single large scale purchase is made. With a better cost trade off between 
t.l.m. components and a more complete view of overall t.Il.m. investment 
compared to other ISS components, a more informed decision can be made as 
to what part of the TSS should be scaled to achieve compatibility between 
them all. 

Concerning the pole error analysis of Chapter 3, several comments 
seem appropriate. It should be noted that none of the pole error analysis 
was carried to more than 39 sections. This is due to the fact that 
none of the root subroutines available was able to generate stable 
solutions beyond this point. Several algebraic root finding methods were 
used in an effort to obtain solutions above the 40th degree, but none 
succeeded. This is believed due to the fact that all the roots are in 
the region between minus one and plus one, and as the degree of the 
polynomial gets higher tine roots tend to crowd near the extremes. Scal- 
ing of the coefficients of the polynomials was attempted but this 
technique resulted in instabilities at even lower degree. Three of the 
four subroutines used in this problem are available in the IBM System/360 
Scientific Subroutine Package (360-CM-03X version III) using three 


different (but algebraic) techniques. The subroutine POLRT uses the 
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Newton-Raphson iterative technique and is explicitly limited to poly- 
nomials of degree thirty-six or less. The subroutine DPRQD uses the 
quotient-difference-displacement method - this routine generally 
became unstable at about the 30th degree for this problem. The DPRBM 
subroutine uses the Bairstow iterative method of quadratic factoriza- 
tion; this routine was the most stable of the four methods - giving 
stable solutions through degree thirty-nine. A fourth subroutine, 
ROOT 1, using Muller's method (written by the computation group at 
M.I.T.) became unstable at about degree twenty-four. The Bairstow 
method is the routine selected for all the programs of Appendix D. 

No extensive effort was made to determine exactly why the subroutines 
became unstable, nor were more elegant matrix methods attempted. It 
was felt that solutions to degree forty gave information sufficient 
for the purposes of this design, so that no further effort was expended 
in this direction. Nevertheless the instability problem 1s at least 


of academic value and worth further study. 
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APPENDIX A 
MODEL TRANSFORMATIONS 


A.|1 Equations of the General Non-Transposed Model 


In the general, non-transposed model of Figure 2.1] Chapter 2, 
the differential equations relating the line to ground voltages and 
the currents are given by Kirchoff's Voltage Law (The Voltage Law): 


aV dl I 


= at (Ry p#Rgg MI, + (Ly y+ gg) FEB + Rogly + (Lyp*Log) ae + 
31 
ane.” (Ly gtlan) a (A-1a) 
~ o>. Peal. beled +13) 2 ile +R..)1, + (LQ,+L,,) “1g 
az 00°a 2 C0 ett 22°00! *b 22°-00/ dt 
31 
Reale (olan) ae (A-1b) 
~ ec. Ri,nl. + (L45+L ) va + Roe Leek 8) b , 
JZ 00°a ee OGee or 00°b 237-00) ot 
31 


(R.+R We + (L,,+L (A-1c) 


aS 
3000 Sons00) ot 
Similarly, the currents can be related to the line to ground voltages by 


Kirchoff's Current Law (The Current Law): 


a. OV av. OV 
az = (Cygtly tly 3) aE - Cyn ae > 813 FE (A-2a) 
al oV ov aV 
Des a ae ae _& 7 
92.7 ~ “2 3 + “aot yatto3) se - 23 3F ae) 
av oV OV 
Coe dy ma) __£& = 
— 9b 7 7 913 at 7 %03 get MO gqtl1 3*C03) ag (A-2c) 


A.2@ Transformation of the General Non-Transposed Model to a “Balanced 
Mutual" flodel 


Equations (A-1) describing Figure 2.1 can be transformed to a 
set describing Figure 2.2 using simple algebra. Let al eal = 3Ip 


and add and subtract 
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31 31 
a b 
Lig at * 443 ae 
3] 3] 
a b 
Lig ot * a3 ae t An 
31 
b Cc 
and Li3 9¢ * 413 3e 
The result is 
3V 31 
a. _ a 
- 97 = Ryytg + (Ly 7443) ae + 
av, al. 
- 577 = (Lyomky3) gE + Rook, + | 
av al. 
- az = (bo3-Ly3) get Rg3l, + C 


Equations (A-2) and (A-3) are the 


mutual" model of Figure 2.2 with: 


bay = bay > 443 Lip 
27 (22 7 73 L53 
L33 = b33 - 443 Li 
Be” “y2 Hy 
Sige ie ie 
O53 * “93 be 


00 ~ 


o3°9 


from eq. (A-1a) 
aL. 
x3¢ ~—sofrom eq. (A-1b) 
from eq. (A-1c). 
a1, 
Liohy3) gE + 3Roglq + 
aly 
3(Logtl 3) a¢- (A-3a) 
a1, a1 
Loon-Ly3) gE + (bo3-Ly3) 5E7 + 
| aly 
BRoglg + Slkogtly3) ge ((A~ 3b) 
a1. 
Lagtli3) a¢~ + 3Roglgg + 
aly 
3(Logtly3) 3¢- (A-3c) 
relations describing the "balanced 
= byo - by; C19 = “0 
= b53 - by3 C50 = “a9 
= log * Ly; C39 = C39 
FS Roo = Roo 
= Roo 
Rag (A-4) 


A.3 Transformation of the General Non-Transposed Model to a Transposed 


Hodel 


One transposition cycle of a transposed three-phase transmission 


line can be visualized physically 


as shown in Figure A.1l. 


In the 
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transposed line let: 


ee = self inductance of a phase = a 
i = mutual inductance between two phases = 
cl 7 : | 
ais Capacitance to ground of a phase = Cy 
= mutual capacitance between two phases = 
R!. = self resi f = RI 
«; = se resistance of a phase = Re 


J 


i 


: 
Cy 


Then from figure (A-1) it is seen that in terms of the parameters of 


Figure 2.1: 
ll = iJ = Ta = rN Meer ee) = J 
11 Ge Sho: 3 S 
J = J = iw = Ele. sme = J 
12 13 23 3 M 
cl = cl = cl = “10 * “20 * “30 = cl 
10 20 30 5 0 
c! = cl = cl = Bie ae es = c! 
12 13 23 3 iM 
A et 1 22 ees a 
11 22 33 3 S 


(A-5a) 


(A-5b) 


(A-5c) 


(A-5d) 


(A-5e) 


The significance of these relations can be seen more clearly by modeling 


a transposition cycle of the transposed line with the model of Figure 


Q> 
aad 
© 


Q OQ 
— ct 


© 


Q OQ 
— ctr 


(A-6a) 


(A-6b) 


2.1. The differential equations corresponding to equations (A-1) with 
I, + I - iP = 31 are: 

: _- Rey + bg =e Li at Ls a Roglo + Shoo 

‘ -- Li = Rol, + bs am Li am = G0g00 > eae 

oe = 1 2 + Ly a Rel, + Lg as oe dag 


(A-6c) 





ie 


Adding and subtracting w Pada a. appropriate in equations A-6 obtain: 
va = Rel, + (Lg - ty a Rog! 5 3 (Ly u Log? “0 (A-7a) 
. ab. Rel, + al - Ly) ab, 3Roglg + 3 (Ly + Log) “0 (A-7b ) 
: Ee - Rte. + (Le = Ly) a + BRgglg + 3(Ly + kgg) - a) 


In the balanced steady-state Ip = Q; the currents in the transposed line 


are positive sequence only, and equations (A-7) reduce still further to 


three decoupled “single-phase equivalent" equations. With Weal 5 = 
positive-sequence inductance and Re = Rp c= positive-sequence resistance: 
OV al. 
- 327 Rp.s.ta * bp.s. oe (A-8a) 
ov. | al. 
- 37 = Rpg ty t bps. 5E- (A~8b) 
aV a] 
C2 eS: Ms 
ES hy MC (A~8¢ ) 
If only zero-sequence currents flow in the line, Is = I, = Tos and 


equations (A-7) reduce to 


av, T i) Aj a 
> ee (URE SI eM ae 2) Sila) ee (A-9a) 
OV 3 
by yall Ht T 0) 
eee Say) ines (Leet cine nD les (A-9b ) 
9V at. 
Coe Ii T T 
as = (Re + BRop)I_ + (Le + 2Ly + 3Lgp) a oe) 


Equations (A-8) and (A-9) describe (in terms of The Voltage Law) the 
symmetrical components equivalent circuits of the transposed line 
described by equations (A-7). The transposed line corresponding to 


equations (A-7) is shown in Figure 2.3. 
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The Current Law equations for the transposed line corresponding 


to equations (A-2) can be seen by inspection of Figures 2.1] and A-1: 


ail QV QV 9V 
aaa T a iis ieee 
32 Xp y 20) aE — Sage > Sn de A= 198) 
a1 aV aV aV 
b _ seh i iin be Real ae 
Same toe oe eo oe Saat Soe 
a1 aV QV 3V 
Gee le ae go T eos 
5207 Sa oe 7 Se ft Mo * 2Cy) oe tans) 
I re a MR ee 
etting V. b - 0 qa acting C,, vs rom 
equations (A-10) as appropriate: 
aI aV 9V 
Tr 
= ape = (Cy + 3Cy) gE - 3G oe oe 
al aV QV 
b i Os T 
= Bz = (Cy + 3Cy) gem = 3G) SEM am 
31 aV QV 
Coe T @ lL es0 
waz 1 0 * 24) oe 7 8Sn Ge cae) 


In tne balanced steady-state Vg = 0; the voltages are positive sequence 


only, and equations (A-11) reduce to three decoupled “single-phase 


ii i 


equivalent’ equations. With (Cp v 3C,,) = Cy Sie positive sequence 


Capacitance, equations (A-11) become: 


ol. OV. 
ag = Sy a (A-12a) 
ol OV 
= Oy : 
ee Gala) 
oT. OVO 
oe pls. GE= (A-12c) 


If only zero-sequence voltages are applied to the transposed line, 


2 a ; oe ij 
UA = Vi, = Vee And with Co 


equations (A-11) become: 


=. C5 5. = Zero sequence capacitance, 
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ol aV 
a _ Saal : 
ince ae (A-13a) 
ol av 
D Yb 
Serpe OS ae (A-13b) 
aq. OV 
Se Gan, AS (A-13¢) 


Equations (A-12) and (A-13) describe (in terms of The Current Law) the 
symmetrical components equivalent circuits of the transposed line 
described by equations (A-11). The transposed line model corresponding 


to equations (A-11) is shown in Figure 2.3. 
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APPENDIX B 
MODEL IMPEDANCE DERIVATIONS 


B.1.1 Derivation of L-Section Model Impedance-Inductor Input- 
Shorted Termination 


Using the m section line shown in Figure 3.3, write difference 


equations near the nth section (0<n<m): 





d] 
Ue - tel = Le dt eer 
qv 
Mreeah ~ I, ~ fe dina 
d] 
7 n+] 
"ntl - ¥, = Le dt i Rone] 


Subtracting (B-1) from (B-3): 


aye, * Va =e 


n-] 


Substitute (B-2) into (B-4): 


dV dV 
Se | 7 Vn u cal 7 LC. dt2 7 Saati : 


Assuming solutions of the form We = Re[V e®*). 


AN 


V 2 


a ZV + ae = § LoC.V, + Ss RCV, 


Next assume V., =K, etn ( 


=A +505) and substitute into equation (B-6): 


et(n-1)A _ 9 inh i: et (ntl )A ae 


-i ) 4 
e -2+e =5 LC. + 5S RoC. 


2 
Sa Ss RC 
Ss Ss Ss Ss 








where X may be a complex constant, 


+nr 
LC. + Ss RC. )e 


(B-4) 


(B-5) 


(B-6) 


(B-7) 





Lo 


Therefore, provided » satisfies (B-7), the solutions to equation (B-5) 


take either of the forms: 


n 


Vn 


v= Re[V eS*] = Re[(K, 


i K aa eo (B-8) 


Re[(A cosh nd + B sinh myer] (B-9) 


To obtain the current, I, use equation (B-1) assuming solutions of the 


form: 
2 7) SIL 
dt ReL I, eo] (B-10) 
= a SIG 
Vi, = RelV, ev]. (B-11) 
Then (B-1) becomes: 
I (s L.+R.) = Viel 
We ; ae 
{ = (B-12) 
n R. +S L. 
In terms of equation (B-9) equation (B-12) can be written as: 
I = A(cosh na - eM e = nda - sinh(n-1)A) (B-13) 


S 


S 


If the short is connected at the zeroth node, then Vo = 0 and from 


equation (B-9) it is seen that A = 0, so that elsewhere on the line: 


Pay 


Vi = B sinh na. (B-14) 
And from equation (B-13): 
* _ B(sinh nd - sinh(n-1)A) 


Finally, the impedance looking into the mtn node - or driving point 
terminals - 1s seen to be: 


(R. + 5s L.)sinh mA 


Z(s) ~ Sinh mA - Sinh(m-1)A 


Where it is noted again that must satisfy 


(B-16) 
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SLC S ROC, 


+ 


cosh A=] + 3 7 (B-7) 








In the case that the model is lossless (Ro = 0), equations (B-16) and 
(B-7) reduce to: 


jul. Sinh ma 


Z (jw) Sinh mA - sinh(m-T)A (B-17) 





cosh (B-18) 


The solution for the poles of (B-17) is indicated in Appendices B.1.2, 
C.2.1, and D.1. 


B.1l.2 Solution for the Poles of the L-Section Model-Inductor Input 
Shorted Termination 


The solution for the poles of equation (B-16) or (B-17) rests 
upon recognition of the fact that the denominator of (B-16) and (B-17) 
1s a polynomial in cosh A of degree (m-1) for all m. That this is 
true can be seen by expansion of the denominator for the first few 
values of m. The expansions of sinh ma up to m = 7 are shown in 
Table B.1. 
Table B.1 


Expansions of sinh ma 


Sinh A = sinh A 

Sinh 2\ = sinh A(2 cosh i) 

sinh 3\ = sinh (4 cosh*¢a- 1) 

Sinh 4X = sinh A(8 cosh*2- 4 cosh i) 


( 
( 
( 
Sinh 5\ = sinh A(16 cosh *A- 12 cosh 
( 
( 


O14] ) 
Sinh 6A = sinh A(32 cosh" i~ 32 och tee COSia) 
sinh 7X = sinhd(64 cosh°a-80 cosh*a+24 cosh@a-1) 





13] 


The expansions listed in the table were obtained by repeated 


application of the identity sinh md=sinh(m-1)\ cosh At+cosh(m-1)A sinh A. 


The use of this identity to obtain sinh md requires the preceeding 


recursion terms sinh(m-1)d and cosh(m-1)A. The term sinh(m-1)d is 


Simply the preceeding term in the table, but cosh(m-1)A must be obtained 


from a separate process of expansions of cosh m\. Table B.2 shows the 


first few expansions of cosh md generated in the loop required to get 


Sinh m,\. 


cosh A 
cosh 2) 
cosh 3) 
cosh 4) 
cosh 5A 
cosh 6A 


cosh 7X 


Table B.2 
Expansions of cosh mA 
cosh A 
2 coshon-1 


Mere ee asin Oh 


aes 


16 Sa See AOE cosh i 


4 2 


8 cosh 1-8 cosh 


= 32 cosh°A-48 cosh A+18 cosh*A-1 


9 


64 ae ROE cosh~A+56 Apso) cosh i 


Using Table B.1, the first seven polynomials set equal to zero in the 


denominator of equations (B-16) and (B-17) are listed in Table B.3. 


Table B.3 


Expansions of Impedance Function Denominator 


1 
] 
] 
= ] 
] 
] 


© © oo. 7] 


cosh A = 0 (m=1 ) 

2 cosh i (m-2) 

2 cosh A-4 cosh*A (m=3) 

4 cosh \-4 cosh¢A+8 cosh? etc. 

4 cosh \-12 cosh2-8 cosh3a+16 cosh’ 

6 cosh A-12 cosh@A+32 cosh3at+l6 cosh*A-32 cosh? A 
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The equations listed in Table B.3 are the equations which must be solved 
to obtain the poles of L-section models of from one to seven sections. 
In obtaining the equations of Table B.3 the sinh A term is divided out 
since sinh A # 0. A more thorough discussion of the constant A is now 


in order. In general, A = d+ ja., so that in equation (B-18), for 











example, 
ee ; WLC. 
cosh X = cosh . cos h. fost de sin hi =]- 5 .  (B-19) 
Equation (B-19) implies: 
wil C, 
cosh A cos A. = 1 - 5 (B-20) 
sinh d Sin h. = 0 (B-21) 
But, again, 
wl C. 
cosh \ = cosh A. cos A. = 1 - 7 (B-23) 


and it is seen that cosh \ is real, and d is pure imaginary (= jA,) if 
2 





W LC. . . . 
a? O. With A = JA,» Ne = 0 so that finally, 
wl C, 
cosh X = cosh jr. = cos h, = 1 - 5 : (B-24) 


Equation (B-24) shows that all values of cosh A must be real and lie 
between minus one and plus one. This, in turn implies that all 
solutions to the "equation" in the denominator of (B-17) must also be 


real and lie between minus one and plus one i.e., 


-1.0 < [cosh A] oe (B-25) 


root — 
But the denominator of equation (B-17) is the same as that of 


equation (B-16), so equation (B-25) applies for the lossy line also. 


( 
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Thus, the equations listed in Table B.3 are polynomial equations in the 
real variable cosh » for which there are numerous "root-finding" sub- 
routines in any good computation facility. With the roots of equations 
like those in Table B.3 readily available, equations (B-7) or (B-16) 


can be solved for s or w respectively to give the poles: 





2 2 
- R / °R 2m 
SCcia, AoE (0) “mae L1-(cosh noes (B-26) 
ana 
Wiss Slee es eee 
ie EE seas (B-27) 
SC L Et root 
Me ers 
Where use has been made of the fact that: 
Me See ame 
Re as Le » and C. = (B-28) 


Equations like those in Table B.3 are seen to be solved easily enough, 
but the problem of how to generate the equations up to fifty or more 
sections (m) remains. The details of the generation problem are left 


to Appendix C. 


B.2 Derivation of L-Section Model Impedance, Capacitor Input, Shorted 
Termination 


Using the m-section line shown in Figure 3.6, write difference 


equations near the nth section (0 < n < m): 


dI 
= n = 
“ntl ~ “n> bs Gt (B-29) 
dI 
ms _ hel - 
et ne] oS. at (B-30) 
dV 5 
Ty Tn os GE era 
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Subtracting (B-30) from (B-29): 


= aed 
Viet 2Yaet “pegs Usama (B-32) 


Substitute (B-31) into (B-32): 





a n 
n 7 Loc. 


n- 1] dt 


Similar to tne development in Appendix B.1.1 assume, 


Vy = Rely, enJety 





n 
7 tnd 
n= Ky e@ 
and obtain: 9 
Ww LC. 
cosh A = ] - 5 (B-33) 
Vo = A cosh nd + B sinh (B-34) 


To obtain the current, Ts assume If = Re[ I. aJuty in equation (B-29): 


jol1, = Vii, - Vy (B-35) 


i sy LA cosh(n+1)A + B sinh(n+1)A - A cosh nA - B sinh nr] 
S 


(B-36) 
With a short at the zeroth node Vo = 0 implies A = 0, and 
V, = B sinh na (B-37) 
2 ae . 
T = Ful, bsinh(nt1)r - sinh na]. (B-38) 


Finally, for the driving point impedance at the mth node, 


: sinh mx z : 
‘sic. ~ Jels Sinn(meT)a-sinh my = 12,3... (B-39) 


and using equation (B-33): 


f 


eS 


Se ll 2 = 
ia fc 0 (cosh ) roots fic = lees ee 


(B-40) 


The solution for the poles of equation (B-39) is indicated in Appendices 


taeee and U.2. 


B.3 Derivation of L-Section Hodel Impedance, Upen Termination 


Using the m section line shown in Figure 3.3, but with the short 


removed, Re = 0, and s = jw, the derivation proceeds as in Appendix 


B.1.1 to obtain: 


ve = A cosn nA +B sinh na 


] 


= +—— Alcosn(n-1)A-cosh nA] + Blsinh(n-1)A-sinh na] 


n Jul, 
But now, I, = 0, and equation (B-42) implies, 


cosh A-] 


B= Sinh A 


Substituting (B-43) into (B-41) and (B-42): 


J =A ar toi a 


n Sinha 
1 = —AL {2 sinh nd-sinh(nt1)A - sinh(n-1)A 
n Jul. Sinh A 


So that looking into the mth section: 


/ 2 jul, F sinh(m+])A-sinh ma 
age 


sinh nt caaditarn iD 
As in Appendix B.1.1, 





wel 
cosh A = ] - = : 
from which 
2 iti 
coc «fl tc Ul - (cosh Nr OGe 


(B-41) 


(B-42) 


(B-43) 


(B-44) 


(B-45) 


(B-46) 


(B-18) 


(B-47) 
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The solution for the poles of equation (B-46) is indicated in Appendices 


teens and 0.3. 


B.4 Derivation of m-Section Model Impedance, Shorted Termination 


Using Figure B.1, write difference equations near the nth 








section: 
dl, 
Vine al s dt. 8-48) 
qv 
I,- 1,20, 4% (B-49) 
dv 
» n+] 
Ik- Ty, 2 Ca . (B-50) 
Substitute (B-50) into (B-48) to eliminate I 
dV 
_ d n+] 
Yn 7 Yast = bs GE Ic, dt” All voi) 
Substitute (B-50) into (B-49) to eliminate I,: 
dV dv 
n+] - n 
L, - . cit baa] 7 C. ‘ie (B-52) 





Write the equivalent to equation (B-51) for the next section toward the 
left: 

d qv, 
-V =L=- ic = + a . (B-53) 


+V], (B-54) 


Substitute (B-54) into (B-51) to eliminate Ge 
Z 

d va 
dt 


dI Fe 
tho ge > bef, a (Vay + V4). (B55) 


y -V 
s dt SS dt 


n n+] 





= LC, 


Rewrite equation (B-53): 


dv, dI 


Ya-1 > Yn = bs°, =” Ls e Brae) 





oy 


L°q aunbey 


(SSa_L$SO]) Lapow UuoLAZdaS-L 
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dT. 
Substitute (B-56) into (B-55.) to eliminate ae 
dv 
Ver 7 2V, + Vig Fac. (B-57) 


Solve equation (B-57) as in earlier cases to obtain: 


2 


cosh A = 1 - w Loc. (B-58) 


and 


V = A cosh nk + B sinh nd. (B-59) 


To find the current use equation (B-56) assuming solutions 


2 7 dwt 
ie = Re{ I. ev] 


Vn 


Re[V, eJuty 


to obtain: 


a 1 . 2 
De TE LV, Bae Ne oe ith LCV J ; (B-60) 
Expanding equation (B-60): 
al] 
mS = aa We Sinh A sinh nd + B cosh nd sinh A]. (B-61 ) 


With a short at the zeroth node, Vo = 0 implies A = 0 so that 


ve = B sinh na (B-62) 
and 
a pg LLY (B-63) 


S cosh nA sinh AX 


At the mth node or driving point terminals : 


e Sinh mA 5 
2, 7 7 Jul, cosh mA sinha ° (B-64 ) 


The solution of equation (B-64) for the poles is discussed in 
Appendices C.2.4 and D.4. The results of the solution are the values 


of [cosh Ae which give the poles of (B-64). Also, using 





SS, 


equation (B-58), 


[1 = (cosh i) 
- root 
cc Cosmin Bae (B-65) 


S=5 


or 
211 - (cosh A aE 


Ste (B-66) 


= 
iN 
25 


S.C. 


Equation (B-66) follows from (B-65) since 


Leeee and Cc. =< 


Som Ss om (B-67) 


B.5 Derivation of m-Section Model Impedance Open Termination 


The derivation of the open circuit case is the same as that of 
the short circuit case through equation (B-61). Now, at the “O"th node, 


lg = 0, and equation (B-61) implies B = 0. Thus 


[= - ar [A sinh A sinh na] (B-68) 
Vv, =A cosh na (B-69 ) 
7 = Vin =a L cosh mA (B-70) 
m I Jue Sinn A sinh mA 


Again, 
- (cosh A 
_m (20) ( Jeet 7 
Yo.c. 2 le (B-71) 


The solution of equations (B-70), (B-71), and (3-1 ) is discussed in 
Appendices C.2.5 and D.5. 
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APPENDIX C 
SOLUTIONS FOR THE POLES OF THE MODEL IMPEDANCE FUNCTIONS 

C.1 Model Impedance Functions 

The model impedance functions obtained in Appendix B for which 
Solutions are sought are shown in Table C.1. For the purposes of 
Appendix C, "solutions" refers to those values of cosh A in the 
denominator of the impedance function which result in the poles of the 
function. In Appendix B.1.2 it was shown that the denominator of 
equation (B-17) is a polynomial in the real variable cosh X. In 
addition it was shown that the polynomial had real coefficients and 


that all solutions for [cosh A] should be such that 


root 


-1.0 < [cosh A] 120 2 


root = 
Since the denominator of equation (B-17) is a polynomial in cosh A, 

it follows that the denominator of all the impedance functions is a 
polynomial in cosh A with all the properties of the polynomial in 

equation (B-17). It is the denominator polynomial set equal to zero 

which results in the poles of the respective impedance function. If 

the denominator polynomials [D(cosh A)] can be generated for any 

number of sections, m, then it is a simple matter to call upon root 
routines in the computer library to solve for the roots of the denominator 
equation, D = 0. An inspection of Table C.] and reference to Appendix 
B.1.2 shows that the "generation" of D involves first the generation of 


Sinh mA and cosh m\ for any m. In Appendix B.1.2 1t was shown that 


sinh mA = sinh A[Polynomial in cosh \] (C-1) 





14] 


Model Driving Pt. Impedance Function 


jul. Sinh mA 


L-Section on ~ §inh mA-sinh(m-1)A (B-17) 
Inductor Input 
Shorted Termination 


{3 


—~ [l-(cosh i) 


be ] (B-27) 


root 


Jul. Sinh ma 


L-Section on = sinh(m+]l) -sinh m_ (B-39) 
Capacitor Input 
Shorted Termination 


= 


wes re L1-(cosh 9 ee (B-40) 


jwL. [sinh(m+1)A-sinh A] 
fn = > (B-46) 
t-Section m 2 sinh mdA-sinh({mt] )A-sinh(m-1 )A 
Open Termination 


we > é~ [1-(cosh ) 


= 


oo (B-47) 


- jul. Sinh ma 
Ln = cosh mA sinh A (B-64 ) 


w= 2 [Fe [1-(cosh A) oe] (B-66) 


~ Jul. cosh ma 


Ly > Sinh A Sinh mA (B- 70) 


m-Section 
Shorted Termination 


m-Section 
Open Termination 


[t-(cosh A) ] (B-71) 


a 
[C 


= 
T 
| 3 


root 


Table of Impedance Functions for Which 
Poles are Sought 


Table C.] 
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and 

cosh mA = [Polynomial in cosh A] . (C-2) 
Denote the polynomial in equation (C-1) by B = B(cosh A) and that in 
equation (C-2) by A = A(cosh i) so that 


Sinh mA = sinh i B (cosh ) (C-3) 


cosh mA = A_(cosh d). (C-4) 


Then it follows that in all the impedance functions of Table C.1 the 
denominator has the form 
D(cosh A) = Linear Combination of A(cosh 4) and B(cosh i). 
(C-5) 
A routine must be found which will enable the computer to generate A 


and Bb combine them to form D = 0, and then solve for the values of 


[cosh lt for use in 
iy ee ondin)>..] (C-6) 
Qo eG root ° 


C.2 Generation of Impedance Function Denominator Polynomial 


In general 
Sinh mA = sinh(m-1)A cosh A + cosh(m-1)A sinh A, 
so that in the terminology of Appendix C.1, 


sinh A B= sinh AB, cosh A +A) Sinh A, 
and 


B,=B =, coshA +A. . (Car) 


Similarly 
cosh mA = cosh(m-1)A cosh X + sinh(m-1)A sinh i, 


so that 
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o ee 
A, = Ap-7 Cosh A + BL. sinh’), 
or 
2 2 
An = Ame] Cosh A + (cosh”A-1)B, 4 . (C-8) 
Throughout the remainder of this appendix denote 
COshsAwoy c- 
and 
Sine byeeese 
so that equations (C-7) and (C-8) become: 
Cs = C Bl if Aad 1) eae ee eee (C-9) 
2 2 : 
av = C Aad +4 C SB, Pa (C-10) 


With m = 1,2...., comparison of equations (C-9) and (C-10) with Tables 
B.1 and B.2 indicates that with Ao = 1 and By = 0 equations (C-9) and 
(C-10) form a recursion loop for generating al] Band A. In 

equations (C-9) and (C-10) it is the coefficients of the powers of cosh i 
which are needed to solve for the roots of the "D" polynomial, so that 

in B and ata it is the coefficients which are of interest. 


In their most general form the polynomials A(m) and B(m) take 


the form (for m = 1,2....): 
a m m- | : 
A(m) = aa mC + a_(m)C ee a, (m) (C-11) 
. m-1 m-2 7 
B(m) = b_(m)C + D1 mde +e. .e. b, (m) ; (C-12) 
Where 
m = number of sections in the model 
C = cosh A 
A(m) = "A" polynomial when model has m sections 


B(m) = "B" polynomial when model has m sections 
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coefficient of cis when model has m sections 


i-] 


coefficient of C when model has m sections 


The subscripts of the coefficients begin at "1" in order to be compatible 


with the subscripting convention in the Fortran IV language. Note from 


Tables B.1 and B.2 that 


1 implies b, (1) = 1, 


i 
— 
w 


C implies a,(1) (C-13) 


and a, (1) = 20. 


To compute B., for example, use equations (C-9), (C-10), and (C-11): 


Bo - 
Where 
and 
In general: 

B(m) 


B(2) = C B(1) + A(1) 
B(2) = C{b,(1)] + Cap(1) C + a,(1)) 
B(2) = [a(1) + b,(1)] ¢ + a,(1) 


= CLb 


+ 


C B(m-1) + A(m-1) 


m-2 m-3 
me] om 1)e +b _o(m-1)C ++. eb, (m-1) J 


[a (m-1)C% +a j(m-1)C™*...e a, (m-1)] . 


Combining coefficients of like powers: 


B(m) 


+ [db y(m-1) + a, (m1)3C" | 


m-2 
+ [b__o(m-1) + ap im-1)C 
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+ Lb (m-1) + a5(m-1) JC 
+ La, (m-1)] (C-14) 


Inspection of equation (C-14) shows that 


b, (m) = Lb; _,(m-1) + a. (m-1) J Zs i S iil (C-15) 

b; (m) = a, (m-1) , i=] (C-16) 
A similar procedure applied for A(m) gives: 

a. (m) = La,_,(m-1) + b, _o(m-1) J i=m, m+] (C-17) 


a.(m) = la;_,(m-1) + Db; _o(m-1)-b, (m-1) J 1=2 a3 eeeeianl 


(C-18) 
a. (m) a La. (m-1) - bo (m-1) J j=2 (C-19) 
a, (m) = [-b, (m-1) J , i=] (C-20) 


Given the initial values for A(1) and B(1), equations (C-15) through 
(C-20) can be implemented in Fortran IV, and the computer made to 
generate and store the coefficients of the "A" and "B" polynomials. 
Having the coefficients of "A" and "B" up to, say, fifty sections 
(m=50), the "D" polynomial can be generated by simple algebra. 
C.2.1 Polynomial for L-Section :lodel, Inductor Input, Shorted 

For the case of the L-Section, inductor input model with a 
shorted termination, the equation for which the roots are sought is 


(from equation B-17) 


Sinh mA - sinh(m-1)\ = 0 


sinh A[B(m) - B(m-1)] = QO. 
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The case X = 0 has no meaning for this model and so 
B(m) - B(m-1) = 0, (C-21) 
and equation (C-21) is seen to be the relation which gives the poles of 
equation (B-17). Thus 
D(m) = B(m) - B(m-1) = 0. (C-22) 
Equation (C-22) must also be written for implementation into Fortran IV 
in the same way that A(m) and B(m) were implemented. Thus: 


ae= d(mc! + ¢ 


il 
Be 
~_— 
= 
— 
] 
Co 
-— 
= 
t 
— 
— 


a np (men ® +....bo(m)C + by (m)] 


m-2 -3 
- [by (m-1)C°* + Bo o(m-1)C" +....b,(m-1) + b)(m-1)]. 
Combining coefficients of like powers: 
a m= | m-2 
D(m) = b (mC + [b_7(m) - bz (m-1) Jc f 


m-3 


[bo (m) - b. ,(m-1)]C ee es 


m-2 
- Lb, (m) - b.(m-1) Jc “- [b, (m) - by (m-1) ]. 


So that identifying coefficients of D(m): 


=: 
~_— 
S 
— 
TH 
o 
— 
= 
— 


- b,(m-1) i = 1]....m-1 (C-23) 
d (m) = b (m) . (C-24) 


Thus the computer can be made to generate the coefficients of the "D" 
polynomial. The form of the coefficients given in equations (C-24) and 
(C-25) are explicitly suited for several "root" subroutines available 
in the IBi-360 Scientific Subroutine Package (SSP) at the Computation 


Center of 11.1.T. Equations (C-15) through (C-20) and (C-23) and (C-24) 
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were implemented in Fortran and solved with the IBIM-360 as shown in 
Appendix D.1]. In Appendix C.1 the roots of D are put into equation 
(B-26) to obtain the poles of a 100 mile long typical 765 kV line. 
The poles of the model thus obtained are compared to the poles of the 
actual line to obtain the pole error plots of Chapter 3. 
C.2.2 Polynomial for L-Section ilodel, Capacitor Input, Shorted 

For the case of the L-section, capacitor anne model with a 
snorted termination, the equation for which the roots are sought is 


(from equations b-39): 


sSinh(m+1)X - sinh m = 0 
sinh ALB(m+1) - B(m)] = 0. 
The case A = UO has no meaning for this model and so 
D(m) = B(m+1) - B(m) = 0. 
Following a procedure identical to that in Appendix C.2.], the 


coefficients of the “U" polynomial are found to be 


d (m) = b,(m) (C-25) 


ae 

“— 
= 

— 
li 


[b,(m)-b,(m-1)]. i=1,2,....m+1 (C-26) 


Note that equations (C-25) and (C-26) appear identical to equations 
(C-24) and (C-23) respectively, but that in the case of equations 
(C-25) and (C-26) the coefficients pertain to a higher power of the 
variable C. Equations (C-15) through (C-20) and equations (C-25) 
and (C-26) were implemented in Fortran and solved with the computer 


as shown in Appendix v.2. 
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C.2.3 Polynomial for L-Section Model, Open Termination 


For the case of the open circuited L-section model, the equations 
for which the roots are sought is (from equation B-46): 
2 sinh m - sinh(m+1])d - sinh(m-1)A = 0 
Sinh mA(1 - cosh A) = 0 
Sinh A B(m)[1-C] = 0. 
Again the case A = 0 adds no solutions, so that division by sinh A can 
be performed assuming X # 0. Finally then 


D(m) = Bi(m)(1-C) = 0 


: -] -2 : 
D(m) = [b,(m)C'” +b, 4 (m)c™™* +....b,(m) + b,(m)](1-C). 

Collecting coefficients of like powers of C: 
D(m) = - b.(m)C" + [b, (m) - b 
+[b.(m) - by (m) JC + b, (m) ‘ 


The coefficients, d., are then seen to be 


d, (m) = b, (m) (C-27) 
d;(m) = [b,(m) - b,_,(m)]  i=2,3,....m (C-28) 
d+] (m) = - b (im) when j=m+] (C-29) 


Equations (C-15) through (C-20) and (C-27) through (C-29) were 
implemented in Fortran and solved with the computer as shown in 
Appendix 0.3. 
C.2.3 Polynomial for m-Section Model, Shorted 

For tne case of the shorted m-section model, the equation for 
which roots are sought is (from equation B-64): 


cosh mA sinh A = 0. 
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Tne case A = Q has no meaning for this model so that 
D(m) = cosh m, 
D(m) = A(m), 
d.(m) = a.(m) | sche ceenaemitl 


This case is particularly simple and is shown implemented in Appendix D.4. 


C.2.4 Polynomial for m-Section lodel, Open 


For the case of the open m-section model, the equation for which 


roots are sought is (from equation B-70): 


Sinh A sinh mA = 0 


sinh?) B(m) = 0 
(1-cosh*) B(m) = 0 
D(m) = (1-C*) B(m) =0 . (C-30) 


Expanding equation (C-30) and collecting like powers of C: 


m+] 


D(m) = b.(m)c™! +B s(m)c™ + [b,_o(m) - b,(m)Jc™! +... 


m- ] 
. Lb, (m) - by (m)]c> “- [b, (m) - b(m) Jc“ - Do(m)C - by (m). 


The coefficients d. are then identified as 


di(m) = -b,(m) if i= 1,2 (cacy 
d-(m) = [b,_4(m)-b,(m)] — i=3,4,...m (C-32) 
d.(m) = b__4(m) j=m+] (a3) 
d.(m) = b,(m) i=m+2 (C-34) 


Equations (C-15) through (C-20) and (C-31) through (C-34) were 
implemented in Fortran and solved with the computer as shown in 


Appendix D.5. 
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